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Abstract. In this paper path integration in two- and three-dimensional spaces 
of constant curvature is discussed: i.e. the flat spaces IR 2 and IR 3 , the two- 
and three-dimensional sphere and the two- and three dimensional pseudosphere. 
The Laplace operator in these spaces admits separation of variables in various 
coordinate systems. In all these coordinate systems the path integral formulation 
will be stated, however in most of them an explicit solution in terms of the 
spectral expansion can be given only on a formal level. What can be stated in 
all cases, are the propagator and the corresponding Green function, respectively, 
depending on the invariant distance which is a coordinate independent quantity. 
This property gives rise to numerous identities connecting the corresponding 
path integral representations and propagators in various coordinate systems with 
each other. 
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1. Introduction 



The invention of the path integral by Feynman [31] is one of the major achievements of 
theoretical physics. In its now 50 years history it has become an indispensable tool in 
field theory, cosmology, molecular physics, condensed matter physics and string theory 
as well [42, 96]. 

Originally developed as a "space-time approach to non-relativistic quantum me- 
chanics" [31] with the famous solution of the harmonic oscillator, it became soon of 
paramount importance in quantum electrodynamics (QED), especially in the develop- 
ment of the nowadays so-called "Feynman rules". It did not take long and Feynman 
succeeded in discussing problems not only in QED but also in the theory of super- 
fluidity. 1 However, it took a considerably long time before it was generally accepted by 
most physicists as a powerful tool to analyse a physical system, giving non-perturbative 
global information, instead of only perturbative, respectively local information, as in an 
operator approach. 

Eventually, a satisfying theory should be based on a field theory formulation, let 
it be the second quantization of the Schrodinger, respectively the Dirac equation, let 
it be a field theory path integral. However, the path integral is quite a formidable 
and difficult functional-analytic object. The very early field theory formulations by 
a path integral, first by Feynman and in the following by e.g. Matthews and Salam 
[35, 82] remained only on a formal level, however with well-described rules to extract 
the relevant information, say, for Feynman diagrams. Indeed, in field theory the path 
integral is cursed by several pathologies which cause people now and then to state that 
"the path integral does not exist". 

What does exists, however, is the very originally Feynman path integral, i.e. 
Feynman's "space-time approach to non-relativistic quantum mechanics". Thanks to 
the work of many mathematicians and physicists as well, the theory of the Feynman 
path integral can be considered as quite comprehensively developed. Actually, the 
theory of the "Wiener-integral" [35] existed some twenty years right before Feynman 
published his ideas, and was developed in the theory of diffusion processes. The Wiener 
integral itself represents an "imaginary time" version of the "real time" Feynman path 
integral. This particular feature of the Wiener integral makes it a not too complicated 
and convenient tool in functional analysis, mostly because convergence properties are 
easily shown. These convergence properties are absent in the Feynman path integral and 
the emerging challenge attracted many mathematicians and mathematical physicists, 
c.f. the references given in [35]. Let us in addition mention Nelson [88] concerning 
the Feynman path integral in cartesian coordinates, DeWitt concerning curvilinear 
coordinates [18], Morette-DeWitt et al. [86] and Albeverio et al. (e.g. [2, 3]) who 
developed a theory of "pseudomeasures" appropriate to the interference of probability 
amplitudes in the path integral. This interference of probabilities, in particular in 
the lattice formulation of the path integral (see below in Section 2) leads to the very 
interpretation of the Feynman path integral. One encounters for finite lattice spacing, 



1 For this kind of general overview we do not want to cite all the relevant references in detail, but would 
instead invite the interested reader to consult some of the textbooks on the Feynman path integral, in 
particular Feynman and Hibbs [32], respectively the collection of Schwinger [101] and the early review 
paper of Gel'fand and Yaglom [35]. 
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i.e. finite N, a complex number $(q l7 . . . , qjv-i) which is a function of the variables q ■ 
defining a q(t), the path integral can be interpreted as a "sum over all paths" or a "sum 
over all histories" 

tf(q",q';T)= £ *[q(f)] = £ e 1 ^/* . (1.1) 

over all paths over all paths 

from q' to q" from q' to q" 

The path integral then gives a prescription how to compute the important quantity 
$ for each path: "The paths contribute equally in magnitude, but the phase of their 
contribution is the classical action (in units of fi). ... That is to say, the contribution 
<&[q(t)] from a given path q(i) is proportional to exp (^<5[q(*)]), where the action is the 
time integral of the classical Lagrangian taking along the path in question" [31]. All 
possible paths enter and interfere which each other in the convolution of the probability 
amplitudes. In fact, the nowhere differentiable paths span the continuum in the set of 
all paths, the differentiable ones are being a set of measure zero (the quantity Aq^ / At,- 
does not exist, whereas (Aq) 2 /At,- does). 

A more comprehensive discussion will be given in a forthcoming publication [58], 
including a table of exactly solvable path integrals, in our lecture notes [59], and 
alternatively, we refer to the several already existing textbooks on path integrals 
Albeverio et at. [2, 3], Dittrich and Reuter [19], Feynman and Hibbs [32], Glimm and 
Jaffe [40], Kleinert [73], Roepstorff [98], Schulman [100], Simon [103], and Wiegel [112]. 
For a short reference we refer to [57]. 

The subject of this paper will be the path integral formulations in (spatial) two 
and three-dimensional spaces of constant curvature. These are the most important 
needed in physics. It is often necessary to consider a given problem from various 
points of views, i.e. coordinate system realizations, say, and moreover find the relevant 
Fourier expansions needed for the harmonic analysis and for the transition from one 
coordinate system to another one, respectively. E.g. in the space A^ 3 ) (the three- 
dimensional pseudosphere) this task was undertaken by Vilenkin and Smorodinsky 
[111] for five coordinate systems on A^ 3 ). Historically, Lame started the project of 
finding all trirectangular systems in Euclidean space which admit separation of variables 
of the corresponding Hamilton-Jacobi and Laplace equations. Further, the works of 
Darboux [16], Eisenhart [26] and Stepanov [105] were set out to solve this problem. 
The result is that there are four real coordinate systems in two dimensions and eleven 
real coordinates systems in three dimensions of the sought type (which e.g. can be found 
in Morse and Feshbach [87]). The corresponding problem for spaces of (non- vanishing) 
constant curvature, i.e. on spheres and pseudospheres was solved by Olevskii [90], and 
a systematic approach for the D-dimensional generalization is due to Kalnins et al. [69, 
70] and references therein. I will return to this classification in the next Section. 

Some of the coordinate systems - cartesian, polar, parabolic - are very familar, 
others much less so. They are all obtained as degenerations of the confocal ellipsoidal 
coordinates (c.f. [84, 87] for such discussions in the three-dimensional Euclidean plane). 
It turns out that the pseudosphere A^ D_1 ^ for a fixed dimension D has the richest 
structure of all of them. This is not too surprising. The uniformization theorem for 
Riemann surfaces states that the fundamental domain of A^ 2 ) /T, where V is a discrete 
fixed-point free subgroups of the automorphisms on A^ 2 ), i.e. V is a Fuchsian group, 
tesselates the entire hyperbolic plane. The Riemannian surfaces may have any genus 
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g > 2 and may be arbitrarily shaped according to the corresponding Teichmiiller space. 
This rich structure makes A^ 2 ) interesting in the Polyakov approach to string theory, 
respectively in 1 + 1-dimensional quantum gravity [42] : in the perturbative expansion of 
the Polyakov path integral one is left with a summation over all topologies of world sheets 
a string can sweep out, and an integral over the moduli space of Riemann surfaces. In 
contrast, the sphere can represent only a genus zero Riemann surface, i.e. a sphere. 

The pseudosphere is on the one hand of particular importance, because the 
manifolds of constant u 2 , u being the four- velocity or the four-momentum, represent 
the physical domain of the variables for particles with real mass (on the mass shell). On 
the other, a similar harmonic analysis as in 

A (2) 

can be studied, c.f. [27, 107], giving a 
tessalation of A^ 3 ^ in the form of three-manifolds. This structure makes A^ 3 ^ interesting 
in 2 + 1-dimensional quantum gravity. 

I am now going to study and compute explicitly as far as possible, all the path 
integral formulations in spaces of constant curvature in two and three dimensions in 
terms of all possible coordinate systems allowing the complete separation of variables. 
Not all our path integral representations will be new. However, we put them all, i) 
into the context of exactly solvable examples of the path integral in a given space of 
constant curvature, ii) We establish in this way numerous identities in terms of a specific 
coordinate formulation of a path integral and its corresponding spectral expansion on 
the one hand, and its explicitly known form of the propagator and the Green function 
in terms of the invariant distance (norm) in the space expressed in these coordinates 
on the other. Spectral theory of functional analysis guarantees the equivalence of the 
identities. In the majority of the cases, especially for where no explicit path 

integral evaluations are possible, this equivalence sets out further progress in finding 
exactly solvable examples for the Feynman path integral by connecting the various 
examples with each other. 

Whereas in this paper we focus on coordinate systems and their corresponding path 
integral representations, one can also look on the matter from a group theoretical point 
of view [10]: For flat space we have the Euclidean group, for the sphere we have 5 , (- D_1 ) = 
SO(D)/SO(£>-l), in particular = SO(4)/SO(3) and = SO(3)/SO(2) and for 
the pseudosphere A^" 1 ) SO (D - 1, 1)/ SO(D - 1). 

In order to look at such a path integral formulation we consider the Lagrangian 
£(x, x) — V(x) (x G M p+q ) as formulated, say, in a not-necessarily positive definite 
space with signature 

(g ab ) =diag ( +!,.. .,+1 , -1,. ..,-! ) ■ (1.2) 

p times q times 

One introduces polar coordinates 

x v = re 1/ (6 1 ,...,6 p+q _ 1 ) , v = 1, . . . ,p + q , (1.3) 

where the e are unit vectors in some suitable chosen (timelike, spacelike or lightlike) set 
[10]. One then expresses the Lagrangian in terms of these polar coordinates and seeks 
for an expansion of the quantity e 2 ( ei ' e2 ) expressed in terms of group elements gi , gi. If 
this is possible one can re-express the path integration of the coordinates x into a path 
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integration over group elements g yielding [10] 



x(t")=x" 

ft 



J X>x(£) exp 

x(t')=x' 



i / ^(x,i)cft ^ / dE x d x J2fLDi n ( 9 ' V) (1-4) 

^ mn fc (1.5) 



Here is defined via the Fourier transformation 

/(<7)= / dE x d x J2fLDi n (g) , / mn = / f^D^ig-^dg , (1.6) 

and (it/ is the invariant group (Haar) measure, f dE x stands for a Lebesque-Stieltjes 
integral to include discrete (f dE x — > as weu as continuous representations. J dE x 
is to be taken over the complete set {A} of class one representations. d x denotes (in the 
compact case) the dimension of the representation and we take 

d x I D^ n (g)DX,* nl (g)dg = 6(\,\')6 m , m ,6 n , n , (1.7) 
Jg 

as a definition for d x . S(X, A') can denote a Kronecker delta, respectively, a 5-function, 
depending whether the quantity A is a discrete or continuous parameter. We have 
furthermore used the group (composition) law 

DLig-'gt) = J2 D k:(9a)DL(9b) ■ (1.8) 

k 

Choosing a basis {b} in the relevant Hilbert space fixes the matrix elements D^ nn 
through D^ nn = (D x (g)b m ,b n ) of the representation D x (g) of the group. In particular 
the L>Q m are called associated spherical harmonics, and the Dq are the zonal harmonics. 
These spherical functions are eigen-functions of the corresponding Laplace-Beltrami 
operator on a, say, homogeneous space, and the entire Hilbert space is spanned by a 
complete set of associated spherical functions -Dom [38, 39, 110]. 

"All what remains" is to look at a convenient representation, i.e. a coordinates 
system. One may say, that this is the very subject of this paper. 

In our case of the free motion things are not too difficult, and the group path integral 
technique can be fruitfully exploited. This kind of approach is always possible in cases 
where a model has a known underlying dynamical group structure (e.g. [65]). Here 
the most famous example is the hydrogen atom with its 0(4) symmetry. Actually, 
this property enabled Duru and Kleinert to apply the so-called Kustaanheimo-Stiefel 
transformation to the path integral problem of the hydrogen atom [23, 24]. This 
dynamical group structure is also important in order to discuss the so-called "super- 
integrable" potentials [30] (where the Coulomb potential and the harmonic oscillator in 
IR 3 are but two examples), and of course their path integral representations; we return 
shortly to this topic in the discussion in Section 6. 
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The further content will be as follows: In the second Section we outline some basic 
information concerning the path integral in curved spaces. We sketch the relevant 
(transformation) techniques along the lines as presented in our earlier work [43, 49, 
55]. Next we sketch the classification of separating variables for Riemannian spaces of 
constant curvature. Here some basic information as found in Kalnins [69] is given. This 
includes the development of separability of variables in the language of path integrals. 

Sections 4 and 5 then deal with the enumeration of the path integral representa- 
tions, spectral expansions, and path integral evaluations in two and three dimensions, 
respectively. These two Sections represent the principal part of the paper and because 
we also list some already known results, some parts will have a review character. 

The last Section contains a summary and a discussion of the results. In the three 
appendices some additional information is given as needed in Sections 4 and 5: In 
Appendix 1 some basic path integrals are cited, in Appendix 2 a particular dispersion 
relation is discussed, and in Appendix 3 the path integral solution on A( D_1 ) in a 
particular coordinate system is given. 



2. Formulation of the Path Integral 

In order to set up our notation we proceed in the canonical way for path integrals 
on curved spaces (DeWitt [18], D'Olivio and Torres [20], Feynman [31], Gervais and 
Jevicki [36], [43, 55], McLaughlin and Schulman [80], Mayes and Dowker [83], Mizrahi 
[85], and Omote [89]). In the following x denotes a D-dimensional cartesian coordinate, 
q a D-dimensional arbitrary coordinate, and x, y, z etc. one-dimensional coordinates. 
A quantity o denotes an operator. We start by considering the classical Lagrangian 
corresponding to the line element ds 2 = g a bdq a dq b of the classical motion in some D- 
dimensional Riemannian space 

Ccfa Q) = y (|) 2 " V(q) = ^ 9ah q a q b - K(q). (2.1) 
The quantum Hamiltonian is constructed by means of the Laplace-Beltrami operator 

■ = -£ A - + - J -L^f^k + v(q) (2 - 2) 

as a definition of the quantum theory on a curved space [94]. Here g = det (g a b) 
and (g ab ) = (gab) -1 - The scalar product for wave- functions on the manifold reads 
[fid] = f dq.^fg /*(q)<7(q), and the momentum operators which are hermitean with 
respect to this scalar product are given by 

V--(— + V -\ V - d J^H (23 ) 

Pa ~ i\dq a + 2 J ' La ~ dq° ' { ] 

In terms of the momentum operators (2.3) we can rewrite H by using the Weyl-ordering 
prescription ([55, 85], W=Wey\): 

H(p, q) = ^(g ab p a p b + 2p a f b p fe + p a p & f b ) + V(q) + AV w {q) . (2.4) 
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Here a well-defined quantum correction appears which is given by [55, 85, 89]: 

AVw = S^ ^^ " R) = 8m" l 9 ^^ + + ^ 



The corresponding Lagrangian path integral reads (MP = Mid- Point): 

q(t")=q" 



K(q",q';T) = ^(q')^(q")] 



v /^'D M pq(t)exp 



l 



£e//(q, q)d* 



q (t')=q' 



W^fq'T^^J^) 1 (II / )1I\ " 



x exp < — 
ft 



m 



^g ab {^)Aq^Aq b - e^q,-) - eAW(q,.) 



2e 



(2.6) 



Here we have used the abbreviations e = (t" — t')/N = T/N, Aq^ = q^ — qj_i, 

q, = 2(Qi + Qj-i) f° r = + J e ) (*j = + ej 5 j = 0, . . . , N) and we interpret the 
limit N — > oo as equivalent to e — > 0, T fixed. The lattice representation can be obtained 
by exploiting the composition law of the time-evolution operator U = exp(— i HT/ft), 
respectively its semi-group property. The Weyl-ordering prescription is the most 
discussed ordering prescription in the literature. 

In an alternative approach the metric tensor is rewritten as a product according to 
g a b = h ac h c b [43]. Then we obtain for the Hamiltonian (2.2) 



h 
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H = -^LB + V(q) = ^h ac p a p 6 h c6 + AVW(q) + V(q) 
and for the path integral (PF - Product-Form) 

*(q",q';T) 

q(*")=q" t „ 
= J y/gV PF q{t)expl^ J 

q(t')=q' 

/ m ■sND/2 r il} r i 



(2.7) 



m . 



KM)h ch (^)q a q b - V (q) - AV(q) 



dt 



N- 



x 



f i N 



m 
27 



/i6c(q J )/iac(q J _i)A^A^ - eV^) - eAV PF (qj) 



(2.8) 



AVpir denotes the well-defined quantum potential 



h 2 

AVp F = — 
8m 



g ab T a T b + 2(g ab T b ), b + g ab >o6 ] + |1 (2/« b % 6 - /* ac , a /i bc , b - ^ c , b /i 6c , a ) 

(2.9) 



arising from the specific lattice formulation for the path integral, respectively the order- 
ing prescription for position and momentum operators in the quantum Hamiltonian. We 
only use the lattice formulation of (2.8) in this paper unless otherwise (and explicitly) 
stated. 

Indispensable tools in path integral techniques are transformation rules. In order 
to avoid cumbersome notation, we restrict ourselves to the one-dimensional case. For 
the general case we refer to DeWitt [18], Fischer, Leschke and Miiller [33], Gervais and 
Jevicki [36], [50, 53, 55, 57], Junker [67], Kleinert [72, 73], Pak and Sokmen [91], Steiner 
[104] and Storchak [106], and references therein. We consider the one-dimensional path 
integral 

X (t")=x" t „ 
K(x",x';T) = [ Vx(t)exp ^[ ( ^x 2 - V(x))dt (2.10) 



x(t')=x' 

and perform the coordinate transformation x = F(q). Implementing this transforma- 
tion, one has to keep all terms of 0(e) in (2.10). Expanding about midpoints, the result 
is 



K(F{q"),F{q')-T)= F\q")F\q') 



-1/2 



m 



lim , 

N^oo \2nien 



N 



1/2 N-l . 

n /<%-n F 'fe) 

j=l J 3=1 



x exp 



h 



N 

E 

i=i 



^F' 2 fe)(A^) 2 - eV{F{q 3 )) - f^^M 
2e vyjA H3> v yH3>> 8m F /4 (^) 



(2.11) 



It is obvious that the path integral representation (2.11) is not completely satisfac- 
tory. Whereas the transformed potential V(F(q)) may have a convenient form when 
expressed in the new coordinate (/, the kinetic term ^F' 2 q 2 is in general nasty. Here 
the so-called "time-transformation" comes into play which leads in combination with 
the "space-transformation" already carried out to general "space-time transformations" 
in path integrals. The time-transformation is implemented [23, 24, 104, 106] by intro- 

rt" 2 

ducing a new "pseudo-time" s" by means of s" = J t , ds/F' A (q(s)). A rigorous lattice 
derivation is far from being trivial and has been discussed by many authors. Recent 
attempts to put it on a sound footing can be found in Refs. [14, 33, 113]. A convenient 
way to derive the corresponding transformation formulas uses the energy dependent 
Green's function G(E) of the kernel K(T) defined by 



- J dTe^ E+ie ^ h K(x",x';T) . (2.12) 



G(x , x ; E) = ( q 



H-E-ie 



where a small positive imaginary part (e > 0) has been added to the energy E. (Usually 
we not explicitly write the i e, but will tacitly assume that the various expressions are 
regularized according to this rule). For the path integral (2.10) one obtains the following 
transformation formulas 



K(x",x';T)= [ ^e- iET /KG(q",q';E) , 
Jtr 2nl 

G(q'\ q'; E) = 1 [F\q")F'{q')] ^ jf° ds"K{q'\ q'; s") , 



(2.13) 
(2.14) 
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with the transformed path integral K{s") given by 

1/2 N-l 



K(q", q'; s") = lim 



m 



n^oo \ 2niefi 



n/ 

7 = 1 J 



f. n r 
J-H ^(A ?J ) 2 -eF' 2 fe)(y(Ffe))-i?) 



3 = 1 



8ml F' 2 (q 3 ) F'{q 3 ) 



(2.15) 



Further refinements are possible and general formulas of practical interest and impor- 
tance can be derived. Let us note that also an explicitly time-dependent "space-time 
transformation" x = F(q,t) can be formulated similarly to the formulae (2.13-2.15), c.f. 
Refs. [53, 58, 59, 73, 93, 106]. 

Finally we consider a pure time transformation in a path integral. We consider 



G(q", q'; E) = y/MW)^ J ds" (q" | exp ( - i a" y/f (H - E) y/f /%) 

(2.16) 

t" 

which corresponds to the introduction of the "pseudo-time" s" = J t , ds//(q(s)) and 
we assume that the Hamiltonian H is product ordered. Then 

• POO ^ 

G{^'^E)= l -{ff") h{1 - D/2) J o *(q",q'; (2.17) 

with the path integral 

q(t")=q" 



Ktf', q'; s") = 



q(t')=q' 



xexp^- 



rn 



-^h ac h ch q a q b - f ( V(q) + AVW(q) - E 



ds } . (2.18) 



Here are h ac = h ac /y/J, \fg = det(h ac ) and (2.18) is of the canonical product form. 

The third ingredient in our calculations will be the technique of separation of 
variables in path integrals [44]. Let us sketch the most important features of this 
technique. We assume that a potential problem V(x) has an exact solution according 
to 



x(t")=x" 



X>x(£) exp 



rn 



x J - V(x) )dt 



x(t')=x' 



= / dE x e- iE - T / h n(x')Vx(x") , 

(2.19) 



and x is a set of variables of dimension d. Now we consider the path integral 

z(t")=z" x(i")=x" 

K(z",z',x",x';T) = | f d (z)G(z)Vz(t) J Vx(t) 



z(t')=z' 



x(t')=x' 



x exp 



m 



(g • z) 2 + / 2 i 2 - 



Z 2 (z) 



+ V(z)+Al/(z) 



eft 



(2.20) 



Here, z denotes a cf-dimensional coordinate with d + d! = D, gi and / the elements 
of the D-dimensional metric tensor g a b = diagf^f, . . . , g#, f 2 , . . . , f 2 ], AV the quantum 
potential of (2.9), and det(g ab ) = f 2d Hi 9 1 = f 2d G{z). As shown in Ref. [44] by 
performing a time-transformation (see Duru and Kleinert [24] and Kleinert [73]) forth 
and back in the path integral (2.20) we can separate the x from the z variables yielding 



K(z",z',x",x';T) = [/(z')/(z 



>/\]-d/2 



dE A ^(x')*A(x") 



X 



|(g.z) 2 -y(z)-Ay(z)--^ y 



dt 



(2.21) 



Of course, also x-depended metric terms can be included (with their corresponding 
quantum potentials) in the separated x-path integration without changing the general 
feature. 



3. Separable Coordinate Systems on Spaces of Constant Curvature 



3.1. The Coordinates Systems 

A systematic approach for the D-dimensional on the problem on the classification of 
separable coordinate systems in spaces of constant curvature (positive, zero, negative) is 
due to Kalnins et al. [69, 70]. By the notion "separable coordinate system" we man any 
coordinate system which separate the classical Hamilton-Jacobi equations, respectively 
the Schrodinger equation. The classification goes at follows [69] : 

3.1.1. The Sphere 5 ,< - D_1 - ) . We denote the coordinates on the sphere g^ -1 -* by the 
vector s = (s , • • • , sd-i)- The basic building blocks of separable coordinates systems 
on g^ -1 ) are the (D-l)-sphere elliptic coordinates 

"eO ' J " " I>? = 1, (3-D 
corresponding to a metric 

D 

{dp,) 2 , P D (p) = H(p-e i ) (3.2) 

(k > curvature). In order to find the possible explicit coordinate systems one must 
pay attention to the requirements that, (i) the metric must be positive definite, (ii) the 
variables {p}^ 1 should vary in such a way that they correspond to a coordinates patch 
which is compact. There is a unique solution to these requirements given by 

eo < Pi < ei < • • • < e D -i < p D < e D . (3.3) 




ds2 ~ 4k^P D (fH 
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3.1.2. The Euclidean Space E^ D \ In D-dimensional Euclidean space we have first the 
coordinate system corresponding to the D-sphere elliptic (3.2) 



■M6IHS' D) (3 ' 4) 



(c 2 constant). In addition there is a second class of coordinate systems, namely the 
parabolic coordinates 

x i = -ji pi H 1" Pr> + e i H H e D -i) , 

„ 2 = c2 rfci(^-^ (3 - 5) 

' II,,/!'' • '.,) 



x =-c ^ t=i x . (j = 2,... 7J D) . 



3.1.3. The Pseudo-Sphere A^ D On the pseudosphere ^ the complexity 

increases considerably. One starts by considering the line element 



D-l 



2=1 ' Jt^ 1 



(dpi)' 



(3.6) 



(/c < curvature), and one must require that ds 2 > 0. It turns out that there are four 
classes of solutions determined by the character of the solutions of the characteristic 
equation Pd(p) = 0. 

3. 1.3. A. The first class is characterised by ^ ej for i,j = 0, ...,D — 1. If 
-D — 1 = n = 2p + 1 is odd then 

• • -Pi-2 < di-2 < Pi-i < ej_i < e; < e i+ i < < e i+2 < • • • < e 2p+ 2 < e 2p+ i , (3.7) 

(i = 0, . . . , p) with the convention that ej, /9j = for j a non-positive integer which give 

p + 1 distinct possibilities. Using = e^+j+i (z = 1, . . . , 2p + 2, j = 1, . . . ,p + 1) and 
e r = e s for r = s mod (n + 1), the coordinates on A^ -1 ) are written in the following 

way 



Ur 



Yl7=i(pi ~ E i 



u, 



Similarly if D — 1 = n = 2p + 2 is even (i = 0, . . . , p) 



(3.8) 



Pi-2 < ei-2 < Pi-i < e;_i < a < e i+ i < pi < ■ ■ ■ < e 2p +i < e 2p 



(3.9) 



3.I.3.B. The second class is characterised by the fact that there can be two complex 
conjugate zeros of Pd(p) = denoted by e± = a + i/5, e 2 = a — i(3 (a, j3 G IR), 
respectively. Together with the convention e^+i = for all other ej there is the one 
possibility 

Pi < fl < P2 < f2 < ■ ■ ■ < Pn-l < fn-l < Pn ■ (3.10) 
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A suitable choice of coordinates is (j = 2, . . . , n) 



(u + iui) 



2 _ i Ui=i(Pi -ct-iP) 



^nr=i(/i-«-^) 



nr=i(pi-/i-i) 



[(«-/i-i) 2 +^]n w -i(/i-/i-i) ' 

(3.11) 



3.I.3.C. In the third class we have the two- fold root e\ = e 2 = a, say. Let us denote 



G 



(0 



0, 



where ej = (7^-2 (j = 3, . . . ,n + 1) with 



gk ^ gi if k ^ I and gk ^ a for any fc, <7 r = p s for r = s mod (n + 1), n = 2p + 1 for 
n odd, and n = 2p for n even, respectively. This case divides into two families with 
coordinates varying in the ranges (i = 0, . . . , p) 

■ ■ ■ Pi-i < 9i-i < Pi < 9i < a < pi+i < g i+ i < ... g n -\ < p n , (3.12a) 
. . . pi_x < gi- t < pi< g t < p i+1 < a < g i+1 < ... g n - X < p n , (3.12b) 

and in either case of n there are p + 1 distinguishable cases to consider. A suitable 
choice of coordinates is 



(u -ux) 2 



= 6 



Ui=i(Pi-a) 
mi\{Gf-a) ' 

9 nr=i(^- a ) 

da mi\(Gf-a) 



(«-Gf2i) 2 n w -i( G « w -^i) 



(0 



(3.13) 



(j = 2, . . . , n) and e = +1 in (3.12a), e = -1 in (3.12b). 

3.I.3.D. The forth case is characterised by e\ = e 2 = es = b. We set ej = hjs 
(j = 4, . . . , n + 1) with hk 7^ &z for 7^ / and hk ^ b for any fc. Then 

• • • < /Oi-l < < Pi < < b < p i+2 < h i+ i < ■■■ < p n -i < h n - 2 < p n , (3.14) 

(i = 0, . . . , p) and there are p+1 distinct cases. A suitable choice of coordinates is 

niLi(Pi-&) 



(u - Ux) 2 
2u 2 (u - ux) 
(u 2 Q -u\- uj) 



u„ 



mzi(Hf-b) ■ 
9 nr = i(p»-&) 

^ 2 nr = i(Pi-&) 
2^ n -j (H w_ 6) 

nr=i(^-^ 2 ) 



(3.15) 



n w - 2 (^ (0 -^ (i -! 2 )(6-^ 2 ) 3 " 

(j = 3, . . . , n) and i/.^ = foj+i (j = 1, . . . , n— 2, z = 0, . . . ,p), and h r = h s mod (n— 2). 
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For the lowest two case we obtain the following small table of distinct coordinate systems 





D = 2 


D = 3 




2 


6 


em 


4 


11 


A (D-1) 


9 


34 



There seems to be no obvious closed recursions for D a natural number [77] . 



3.2. Separability of Variables in the Path Integral 

Separation of variables is always of great importance in many calculations. The 
separation formula (2.21) provides us already with a prescription in the case where 
we know a separating coordinate system. However, the situation is often somewhat 
more complicated and one has to look first for the coordinate systems which separate 
the relevant partial differential equations, i.e. the Hamiltonian, and, more important 
from our point of view, the path integral. In order to develop such a theory we consider 
according to [87] the Lagrangian £ = ^ "Y^h=i ^1^1 an d the Laplacian Alb, respectively, 
in the following way (where only orthogonal coordinate systems are taken into account) 

Alb = j2 1 9 fnJLiMtf}) d 



I 
D 



• X^TTd — TTTTTT^F \9i(€i, • • • ,&-i,6+i, • • • >£d)/(&)-^t~ ) , 

fenf=ify v ^ (3.i6) 

where {£} denotes the set of variables (£i, . . . , £r>), and the existence of the functions 
fi,gi is necessary for the separation [69, 87]. Note that the factorizing in terms of the 
hi is perfectly adopted to the product form prescription. We introduce the Stackel- 
determinant [69, 87, 90] 



D 



MUD = det($y) = H , Mj^i,...,^,^!,...,^) = = h 2^^ , 

(3.17) 

and abbreviate T; = Then 



D 



= M i (£ 1 ,...,&_ 1 ,& +1 ,...,£ D ) n fj&) . ( 3 - 18 ) 

which fixes the functions <7j. Introducing the (new) momentum operators 
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we write the Legendre transformed Hamiltonian [54] as follows 



h 2 

= H - E = A LB - E 

2m 

-2 d 



1 d (Uk=i k k d 



2m 



i=l 



2m 5 



! r; ( 



fidb vac 



E 



ti 2 l 



2m S * 



5> 



1 

S 



i=i 

D 



d 2 



d 



dt-i l dC 
h 



-E 



2 D 



2- E M^- B5 + |- E M,(r? + 2r;) 

i=l i=l 



(3.20) 



Therefore we obtain according to the general theory the following identity in the path 
integral by means of the space-time transformation technique 



D 

n 

i=l 



hiD£i(t) exp 



n 



m 
~2 



D 



J2hHf-^v PF ({0) 



i=l 



dt 



i=i 



S_ 



V£i(t) exp 



ft 7*' 



m 
2~ 



(5'5")^ (1 - fl/2) / 
Jtr 



dE c -iET/h 



i=i 



2yrft 



/■OO ^ 

/ *"n / M r 1/2 ^ 



x exp 



2 D 



f£|+^-|;X>(r ? + 2r< 

i=l 1 i=l 



(is 



(3.21) 



4. Separation of Variables in Two Dimensions 

The subject of the next two Sections will be the path integral formulation in the 
three two-dimensional spaces of constant curvature, the flat space E^ 2 \ the sphere 
and the pseudosphere A^ 2 \ respectively the path integral formulation in the three 
three-dimensional spaces of constant curvature, the flat space E^ s \ the sphere S^ 3 \ 
and the pseudosphere A^ 3 ). For notation, in the two-dimensional case we denote 
by x = (x,y) = (xi,X2) the flat space coordinates, by s = (so,si,S2) coordinates 
on the sphere, and by u = (uo, u±, U2) coordinates on the pseudosphere, and in 
the three-dimensional case x = (x,y, z) the flat space coordinates and s and u in 
obvious generalization from the two-dimensional case. In each case, we first state 
the general form of the propagator and the Green function in D = 2 and D = 3, 
respectively, in terms of the invariant distance (norm) d E (n) (q", q') = x' • x" in E^ D \ 
cos^ S ( C -i)(q // ,q / ) = s' • s"/R 2 in S^" 1 ), and coshd A(C -i)(q", q') = u 7 • u"/R 2 in 
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^(D-i)_ gp^gj-g g(D-i) ig deggribed by the constraint s 2 = R 2 , and A^ -1 ) by 

u 2 = Uq — XliLi 1 u i = R 2 which fixes the signature of the metric in u-space. The metric 
is determined via the classical Lagrangians, say, 

rgrw) = , cgr\s,s) = ^ , ^-"ku) = -f . <4.d 

We put R 2 = 1, since is just a scaling factor of the systems. Consequently, the 
Laplacians are determined by 

D d 2 0-1 d 2 d 2 D_1 d 2 

= E ' a s(-d = E > A A( r.-i) = ^-E^j • (4 - 2) 

1=1 ' i=0 ' u i=l 1 

We would like to point out that the various path integral evaluations are straightforward 
by applying the basic path integral solutions of Appendix 1. For completeness we present 
in this Section some examples somewhat more comprehensive to help understanding; in 
the next Section it is then sufficient just to take reference to the former ones. It is obvious 
that in two dimensions, the more easy path integral solutions are already known and the 
emphasize of this Section will be more on the collection of the various identities, instead 
on the path integral evaluations. On the other, where the theory of special functions is 
not well developed, such path integral evaluations cannot be explicitly done, therefore 
only an indirect reasoning is possible. Especially in these cases we state the general 
form of the propagator, respectively the Green function. We do not discus coordinate 
systems which are equivalent with a presented one and which can be obtained by a 
simple reparameterization of the variables. Compare e.g. the two different definitions of 
the parabolic coordinates [30, 87] in E^ 2 \E^\ and of the equidistant coordinates [45, 
90, 111] for A^ 2 ), respectively. 



4.1. The Flat Space 



4.1.1. General Form of the Propagator and the Green Function. First of all we discuss 
the general form of the propagator and the Green function in the two-dimensional flat 
space E( 2 \ which are easiest calculated in cartesian coordinates (see below). One finds 
that in E^ they are four different coordinate systems [87]. In terms of the norm 
d E (2)(q", q') (where q denotes any coordinate system) they are given by 

K E{2) (d E(2) ( q ", q ');T) 
G Ei2) (d E M'^)-E) 

The norm in E^ is given by (for the definition of the coordinates see below) 
Cartesian Coordinates: 

4 (2) (q",q') = |x"-x'| 2 , (4.5a) 



m 



lm 



2nihT 

m 



6XP 27iT C ^ (2) ( q " ,q/ ' 1 



(4.3) 
(4.4) 
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Cylindrical Coordinates: 

= r' 2 + r" 2 - 2r'r" cos(</>" - 0') , 
Parabolic Coordinates: 



1 r 



(v" 2 + r 2 ) 2 + (v' 2 + e 2 ) 2 - 2(v /2 - e 2 )W' 2 - r z ) - 8 V va 



Elliptic Coordinates: 

= ^(cosh^' cosh \i" cos v' cos v" + sinh^' sinh//' sin v' sin v") . 



(4.5b) 

(4.5c) 
(4.5d) 



4.1.2. Cartesian Coordinates. We consider the usual cartesian coordinates (x, y) = 
x G IR 2 . The metric is given by (g a b) = I2 and for the momentum operators we have 
p x = — i hd x , p y = — i hd y . Therefore: 



h 2 A h 2 ( d 2 d 2 \ 1 , 2 2 . 



2m ^ ' 2m\dx 2 dy 2 ) 2m 
The path integral formulation is well-known and given by [31, 32] 

x(t")=x" 



(4.6) 



x(t')=x' 



PxWexpl-^ ^j = 2lifiT ~, V2fir 



exp(; 



x — X 



IR 



dp 

2 4tt 2 



exp 



iftT 
~2m 



p^ + ip-(x"-x') 



(4.7) 



(4.8) 



Note the formulation via a path integration over the Euclidean group [11]. 

4.1.3. Cylindrical Coordinates. We consider two-dimensional polar coordinates 

x = r cos 4> > r > , 

y = r sin , < < 2/T . 



(4.9) 



The metric is given by (g a b) = diag(l, r 2 ), and the momentum operators have the form 



h ( d 1 



i V dr ^ 2r 



This gives for the Hamiltonian 

n 2 A / a 2 1 a 1 a 2 

-A B (2) = -— ( — + - — 



f>0 



(4.10) 



1 f 2 1 2^ V 



2m &v ' 2m\dr 2 r dr r 2 dcp 2 J 2m\^ r r ^ rq> J 8mr 2 
We therefore obtain the path integral identity [4, 25, 54, 92, 104] 

r(t")=r" 4>{t")=4>" 



(4.11) 



/ rVr(t) J ^(OexpUjf |(r 2 + r 2 </> 2 ) + 



r(t')=r' 



Ht')=<t>' 



8mr 2 



dt 
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r(t")=r' 

Ai(<t>"-<t>') r \ /•*" / 

= M-^EV / ^ (i)exp U ( 



r(t')=r' 



r2 _ 1 

m ;2 *2 1 4 



-r 2 - /V 
2 2mr 2 



•,U" rh'\ r{t " )=r " 
e ll(<t> -<t> ) 



(rVr l/ 2 £ 



ze^ 



2tt 



r 2 ] exp 



i m 



2Ti J H 



r 2 dt 



r(t')=r' 



m 


exp 


im 


2nihT 


_2hT 


m 


exp 


im 


2k\KT 


2hT 



I i //2\ 

-(r + r ) 



" t I mr'r" 



(> 2 + r " 2 _ 2r 'r" cos(</>" - 0')) 



ze^ 



2^ 



/•OO 

/ pdp Ji (pr')Ji (pr") e~ 1 
Jo 



hp 2 T/2m 



(4.12) 

(4.13) 
(4.14) 



4.1.4. Elliptic Coordinates. We consider the coordinate system 

x = d cosh fx cos z/ , fi > , 

y = d sinh sin z/ , — 7r < z/ < 7r . 



(4.15) 



The metric is (<7 a &) = <i 2 (sinh 2 jU-f-sin 2 i^) H 2 and we obtain for the momentum operators 

h ( d sinh fj, cosh \i \ h ( d sin v cos v \ 

Pv = -[^- + —T2 ~ 5— , Pv = -\TT + —T2 ~ 5— • ( 4 -16) 

1 \o/j> sinh + sin uj \ \ov smh /i + sin uj 

Consequently we have for the Hamiltonian 

_h^ A = h 2 f^L + ^_ 

2m E(2) 2md 2 (sinh 2 \x + sin 2 v) \d\x 2 dv 2 

= ^-*-^(sinh 2 fj, + sin 2 z/)~ 1/2 (p 2 + p 2 )(sinh 2 /i + sin 2 z/)" 1/2 



(4.17) 



The path integral construction is straightforward, however no explicit path integration 
is possible. Actually, an expansion into the corresponding wave-functions in the 
coordinates fi and v yields the Mathieu functions me^(?7, h 2 ) and Me^((,/i 2 ) (h 2 = 
mEd 2 /2ft 2 ), respectively, as eigen-function of the Hamiltonian, a specific class of higher 
transcendental functions [84]. However, because we know on the one side the eigen- 
functions of the Hamiltonian in terms of these functions [84], and on the other the 
kernel in E^ in terms of the invariant distance d E {2) we can state the following path 
integral identity (note the implemented time-transformation) 

J Vfi(t) J Vv(t)d 2 (sinh 2 /j + sin 2 v) exp T^-d 2 J (sinh 2 fj, + sin 2 z/)(/i 2 + v 2 )dt 

H(t')=n' v(t')=v' '- 
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m(0=m" 



'IR 



dE 
2nK 



v(s")=v" 



i ET/h 



J ds" J V/i(s) J Vu(s) 



xexp^- 



m 



^(/i 2 + z> 2 ) + Ed 2 (si-nh z \i + sin 2 v) 



v{0)=v' 

■2/„;„l_2 



ds 



(4.18) 



- V 

2tt ^ 



pdp e" 



ihp 2 T/2m 



X 



^) me,(V, Me« *(£', ^) Me?^", ^) 



m 
2nihT 



exp 



lm 
2ftT 



4( 2 >(q",q') 



(4.19) 
(4.20) 



and <i E (2) (q", q') must be taken in elliptic coordinates. The functions me^ry, h 2 ) and 
Mej, 1 '^, /i 2 ), are mutually determined through the separation parameter A = \ v {h 2 ) 
yielding a countable set of numbers v e A. In particular, the functions Me^(z, /i 2 ) 
yield in the limit /i 2 — > the Bessel functions J^, i.e. Me^^/Zi, /i 2 ) ~ ^(2) — > 0), 
and the functions me^z, /i 2 ) oc e luz (h 2 — > 0), therefore obeying the correct boundary- 
conditions of our problem [84] . 



4.1.5. Parabolic Coordinates. We consider the coordinate system 
x = tV , 2/ = ^ 2 "£ 2 ) , ^GlR,r ? >0 



(4.21) 



(alternatively £ > 0, n e IR [84]). We have (<7 a &) = (£ 2 + ry 2 ) ]1 2 , and consequently for 
the momentum operators 



n ( d 



i £ 2 + ?7 2 

This gives for the Hamiltonian 



h ( d 



Pv = T [-ZZ + 



rj 



i \drj £ 2 + rf 



(4.22) 



2m 



A E (2) = 



(&_ 

2m(£ 2 + ?7 2 ) \d£ 2 



d 2 
drj' 



(4.23) 

Note AVpp = which is a peculiarity of two dimensions with metric oc 1. These 
coordinates have been used in literature to discuss the two-dimensional "Coulomb- 
problem", c.f. [24, 52, 64]. The transformation (x,y) i— > (£,rj) actually is the two- 
dimensional realization of the Kustaanheimo-Stiefel transformation. The arising path 
integral in these coordinates must be treated by a time-transformation due to the 
metric terms in the Lagrangian. The transformation in its continuous and lattice 
implementation, respectively, has the form 

" w = i pwtW) • ' = ( ^» As ' • < 424) 

and decouples the £- and ry-path integration giving two harmonic oscillator path integrals 
in £ and r\ with frequency to = yJ—2E/m, respectively. Because E > in our case, 



17 



u> is purely imaginary and we obtain two repelling oscillators. With the help of the 
harmonic oscillator path integral and the relations in Appendix 2 we obtain the path 
integral identity 



*(*")=€" 



v(t")=ri" 



I ^ I 

€(*')=€' v(t')=v' 

dE 



X>7/(*)(r + ^)exp 

««")=*" *?(»")=»/' 



1 m 
2ft /„ 



(e 2 + ^)(^ 2 + r) 2 )^ 



IR 



2yrft 



t-iET/h 



V V (t) 



€(0)=€' 



»7(0)=77' 



xexp^- 



y(^ 2 + r) 2 )+^ 2 + ^ 2 ) 



(is 



x exp 



e- iET ^ / <te 



moo 



tv i ft sinews" 



1 moj 



2ft sincus 



77 +£ +^7 +^ Jcosws 



cosh 



mo)(^" + f/V n 



iftsincus" 



(4.25) 
(4.26) 



and denotes the summation over even and odd states respectively; the functions 
(S^) are given by 



%,C «»W ^2 



r/2ap 



|r( 



1 KM 2 p(0) 



X 



4 2P /| ^| +iC/p (e- i7r/4 v^0^_ iC/p (e- i7r / 4 v^^) 



r (f-i)| 2 ^ + i C/P ( e " i7r/4 v / ^e)^ 1 |_ iC/p (e- i7r / 4 v^r 7 ); ' (427) 



which are ^-normalized according to [79] 

[°°d V [ rfe(e 2 +^ 2 )*^*(e^)<T ) (e^) = ^ / -p)^(c / -c) • 



(4.28) 



Note that in the evaluation of the path integral one has to take into account that by 
using the harmonic oscillator solution for the £- and r\- variable, respectively, one actually 
uses a double covering of the original (x, y) = (xi, x%) G IR 2 -plane, i.e. u = (£, 77) e IR 2 . 
Furthermore we have taken into account that our mapping is of the "square-root" type 
which gives rise to a sign ambiguity. "Thus, if one considers all paths in the complex 
z = x + iy-plane from z' to z" , they will be mapped into two different classes of paths in 
the u-plane: Those which go from u' to u" and those going from u' to — u". In the cut 
complex z-plane for the function |u| = -\f\z\ these are the paths passing an even or odd 
number of times through the square root from \z\ = and \z\ = —00. We may choose 
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the u' corresponding to the initial z' to lie on the first sheet (i.e. in the right half u- 
plane) . The final u" can be in the right as well as the left half-plane and all paths on the 
z-plane go over into paths from u' to u" and those from u' to — u" " [24]. Thus the two 
contributions arise in (4.25). The last line of (4.26) is then best obtained by considering 
the Coulomb problem —qxqijr in IR 2 , applying (A. 2. 2), respectively (A. 2. 4), performing 
a momentum variable transformation (p£,p v ) — > (^(a + 0) ~~ 0) ( a = ^ 2 / m( ?i ( ?2 
is the Bohr radius) with the new variables (p, () and finally setting (71(72 = 0, i.e. a = 00. 



4.2. The Sphere 



4.2.1. General Form of the Propagator and the Green Function. The sphere is 
the first non-trivial space of constant positive curvature. There exist two coordinate 
systems which admit separation of variables on S^: polar coordinates and sphero- 
conical coordinates. The invariant distance cosifj S (2) = s' • s" in these coordinates is 
given by (for the definition of the coordinates see below) 



Spherical Coordinates: 

cos ip s( 2) (q", q ) = cos 9" cos 9' + sin 9" sin 9' cos(0" - <p') , (4.29a) 

Sphero- Conical Coordinates: 

= sn n" sn [i dn v" dn v' + cn [i' cn \i cn v" cn v' + dn \i' dn \i sn v" sn v' . 

(4.29b) 

The propagator and the Green function on are best calculated in terms of polar 
coordinates. One obtains 



K s(2) (s^ (q", q'); T)=J2 ( cos ^ (3 , (q", q')) exp 



G sW U sm (<i"^);E 



1=0 



m 



ihT 

~2m 



1(1 + 1) 



i + V ^ 2 +i/4 ( " cosV5( 2 )(q", q')) 



2h' 



sin 



n(\-^2mE/h 2 + \) 



(4.30) 
(4.31) 



4.2.2. Spherical Coordinates. We consider the polar coordinates 

so = sin 9 cos <fi , < 9 < tt , ) 

si = cos# , < (/>< 2n , > (4.32) 

S2 = sin6>sin</> . J 

These are the usual two-dimensional polar coordinates on the sphere. The metric tensor 
is (g a b) = diag(l, sin 2 9), and the momentum operators have the form 

% ( d 1 \ h d , inn . 

p ° = l(Te + 2 cote ) ■ P * = T^' (433) 
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For the Hamiltonian we obtain 

h 2 , 



2m 



h 2 ( d 2 B Id 2 

cot 0— + ■ 



2m V d0 2 



1 



86 



sin' 



i 



2m\ Pe + sm 2 l '"J 8m ^ ' sin 2 
The corresponding path integral is well-known [10, 55, 92] and we have the identity 

9{t")=9" 4,(t")=<t>" 

J sm0V0(t) J V<f>{t) 
9{t')=9' 4>(t')=4>' 



^11+ 



1 



(4.34) 



XeXP H./, 



m 



^-(e 2 + r 2 sm 2 e<j ) 2 ) + ^-(i + - 1 



8m V sin 2 6 
P f> + l)Pi ( cos ^) (q", q)) e" 1 ^ (Z+1)/2 



dt 



47T 

oo Z 



= J2J2 yr*(o',<i>')Yr(o",<i>")e~ ihTl(l+1)/2rn ■ 

1=0 m=-l 

The Yl n {6,phi) are the usual spherical harmonics on the S^-sphere. 



(4.35) 
(4.36) 



4.2.3. Sphero-Conical Coordinates. We consider the general elliptic coordinate system 

(4.37) 



4 + s\ + s 2 2 = 1 



S , S l , S 2 



+ 



+ 



Pi - a pi - b pi - c 



, 



where c < p2 < b < pi < a. The connection to the variables s on the sphere is 
given by 



2 _ (pi -a)(p 2 -a) „ 2 _ (Pi ~ & )(P2 - 6) „ 2 _ (Pi ~ C )(P2 - c) 

(c — a)(c — 6) 



(a-c)(a-6) Sl (b-a)(b-c) 



> s 2 



(4.38) 



One now can make the identification pi = /c 2 cn 2 (p, k), p2 = — /c /2 cn 2 (z/, &;'), 6 = — 
a = k 2 and c = 0. This yields 



s = sn(p, fc) dn(z/, A;') , k 2 + k' 2 = 1 , 

51 = dn(p, k) sn(z/, /c') , 

52 = cn(p, k) cn(^, fc') . 



(4.39) 



The functions sn(p, k), cn(p, k) and dn(p, k) denote Jacobi elliptic functions, where e.g. 
cn(p, k) has the periods 4 i K, 2K + 2 i K', 4K, respectively (K, K' the elliptic integrals) 
[1]. Note the relations cn 2 p + sn 2 p = 1 and dn 2 p = 1 — k 2 sn 2 p. The metric tensor 
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g a b in these coordinates is therefore given by (g a b) = (/c 2 cn 2 / u + k' 2 cn 2 v) 1 2 . The 
momentum operators are 



h ( d k 2 sn /j, cn fx dn 
Pm i V^A* /c 2 cn 2 (j + k' 2 cn 2 z/ 



h ( d k' sn v cn z/ dn z/ . 
i k 2 cn 2 + fc' cn 2 v 



and the Hamiltonian has the form 



2m 



A S (2) = 



d 2 d 2 
+ 



2m k 2 cn 2 /i + k' 2 cn 2 v \ dfi 2 dv 2 



(pI+pI: 



^ m a/ k 2 cn 2 \i + k' 2 cn 2 v 



yj k 2 cn 2 \x + k' 2 cn 2 z/ (4.41) 



The path integral can therefore be written down yielding [50] 



/*(*")=*»" 



**(*')=/*' 



y y Vu(t)(k 2 cn 2 n + k' 2 cn 2 u) 

t" 

[ (k 2 cn 2 Li + k' 2 cn 2 v)(/l 2 + v 2 )dt 
Jv 



x exp 



lm 

2ft 



K 1=0 

oo 



= i E( 2Z + W ( cos ^ (q", q')) e" iftT,(,+1 



)/2m 



(4.42) 
(4.43) 



z=o 



The eigen-functions ^a(a*) ^) = Ai jR (/j,)Bi jK (v), are Lame polynomials so that the 
functions Ai jK and _E?z )K are solutions of the Lame differential equations [28] 



d 2 A; ;K (/x) 

d^ 2 

d 2 £?z, K (z/) 
dz/ 2 



+ 



+ 



-^ + /(/ + l)-/(/ + l)A; 2 sn 2 ^ 



AM = , 



j - 1(1 + l)k' 2 sn 2 v 



(4.44) 



Bi, K [y) = . 



These functions are highly transcendental ones and cannot be expressed in terms of 
standard power series expansions like (confluent) hypergeometric functions. 
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4.3. The Pseudosphere 



4.3.1. General Form of the Propagator and the Green Function. The evaluation of the 
Green function on the Pseudosphere A^ 2 ) was subject of a couple of papers. It was first 
evaluated by means of path integration in [54], followed by [45, 56] (compare also [60, 
76]). The propagator and the Green function have the form 



K A(2) (d AW (q", q'); T) = J pdp tanhyrp 
xP ip _i (cosh<i A (2)(q // , q')) exp 
G A (d A ( 2 )(q",q');£) = -^Q_ 1/2 _^ 2rnE/h 2_±_ (coshrf A(2 )(q // ,q / )) 



i HT ( 2 1 



(4.45) 



(4.46) 



The invariant hyperbolic distance in A^ 2 ) is given by (for the definition of the coordinates 
see below) 

Pseudo-Polar Coordinates: 

cosh d A (2) (q", q') = cosh r" cosh r' — sinh r" sinh r cos(0" — <j)') , (4.47a) 
Horicyclic Coordinates: 

(*"-*') 2 + y"W (4 47b) 

Equidistant Coordinates: 

= cosh^' — r^) cosh r" cosh r{ — sinh r" sinh r{ , (4.47c) 

Pseudo-Ellipsoidal Coordinates: 

= nc p" nc p' nc i/' nc v' — dc dc /i' sc i/" sc v' — sc sc p' dc i/' dc z/ . 

(4.47d) 

Because the discussion of the path integral solutions on the pseudosphere has been 
done in some extend in Refs. [45, 54], we just cite the results, with the exception of the 
case of the ellipsoidal coordinates. In the nomenclature we follow [69] . 

4.3.2. Horicyclic Coordinates (Poincare Upper Half-Plane). We consider the coordinate 
system 

u o = 7r{y 2 + 1 - x) , U! = - , y > , 

2y y 

u 2 = —(y 2 - l - x) , x e m , 

^y 

The metric has the form (g a b) = and the momentum operators are 

fc = TS- p » = lU-?j' (449) 

22 



(4.48) 



For the Hamiltonian we have 

-2 t 2 



2m 



h 2 2 { d 2 d 2 \ i, 2 2 2 , , Arn . 
A ^ = 2^ y \^ + W 2 ) = ^ iy pl + ypiy) • (4 - 50) 



We therefore obtain the path integral identity [45, 54, 60] 

x(t")=x" y(t")=y' 



I Vz{t) I 



x(t')=x' 

vW 7 



y(t')=y' 

f oo 



z>y(*) 

— o— exp 



/»/ «9 • 9 

lm / ar + y z 



dt 



2h J t , 

^f- [ dk [ pdp smhirpe-' lhTip2 + ^ /2m e ikix ''- x,) K ip (\k\y')K ip (\k\ 

71 JlR Jo 



(4.51) 



4.3.3. Equidistant Coordinates (Hyperbolic Strip). We consider the coordinate system 
uq = coshri coshr2 , u\ = sinhri , i~i,T2 G H , 

«2 = cosh 7~i sinh T2 . 

The metric has the form (g a b) = diag(l, cosh 2 ri), and the momentum operators are 



(4.52) 



nf d 1 



Pn = - 



T1 i V^n ' 2 

For the Hamiltonian we have 



+ - tanh T\ 



h d 



h 2 ( d 



2m 



d 



A A(2 ) = - — h tanh n- h 



1 <9r 2 

1 a 2 



(4.53) 



2m \ <9ti 
= 2^ 1 P - + 



cosh 



<9ti cosh n 9r 2 

1 2 \ L 

ti 2 r/ T2 J + 8mV cosh 5 



Tl 



(4.54) 



Therefore we obtain the path integral identity [45] 



(ilt\ II It 

Tl (t )=T 1 r 2 (t )=T 2 



coshriPri(t) 



Vr 2 {t) 



ii f m ( -i ,2 -2\ ft 

( COsh T[ COsh r ^" 1/2 e-i^/ 8T n I g e i/e(T»-r 2 



1 



cosh Ti 



Tl(t )=T 1 



X>ri (t) exp 



n(t')=Tj 



m . 2 h 2 k 2 + \ \ 



hj t , V 2 Tl 2m cosh 2 n 7 
( cosh r( cosh rfr V2 / ^e 1 ^'"^ 



1 



p sinh 717? cip 



2 Jir cosh 7r/c + sinh 717? 



e - i^T( P 2 + i)/2m p - ip ( tanh r {) J R i P_ i (tanh r{') 



(4.55) 

In the path integral evaluation one successively uses the path integral solution of the 
special case of the modified Poschl- Teller potential (the simple Rosen-Morse potential 
[46]). 
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4.3.4. Pseudospherical Coordinates (Upper Sheet of the Two-Sheeted Hyperboloid, 
Poincare Disc). We consider the coordinate system 



uq = coshr , t > , 

u\ = sinh r sin </> , < </> < 2/T , 

U2 = sinh r cos 4> . 



(4.56) 



These are the well-known two-dimensional (pseudo-) spherical coordinates on A^ 2 ). The 
metric has the form (g a b) = diag(l, sinh 2 n), and the momentum operators are 



hf d 



Pr = 



i\dr ' 2 
For the Hamiltonian we have 



H — coth r 



V<\> 



h d 



h 2 fd 2 



d 



2m 



A A(2 ) = - — — + cothr — + 



i d<t> 
1 d 2 



(4.57) 



2m \ dr 



dr cosh r dcj) 



2m 



Pr + 



sinh r 



P<t> + 



8m 



1 



sinh Ti 



(4.58) 



Therefore we obtain the path integral identity [45, 56] 

t(*")=t" 4>{t")=4>" 



I 



smh.TVr(t) 



Vcf>(t) 



T(t')=T' 



x exp 



i 



m 



(f + sinh r<; 



8m 



1 - 



1 



sinh r 



= — ?r / dp y p sinh 7rp exp 

2 - 2 ^ 



ihT 

~2m 



p2 + 4 



F[-+ip + l 



e i l(4>"-(f>') v -i i (cQsh r ')7>-i 1 ( cosh r ") 



(4.59) 



(Pseudospherical polar coordinate system on the disc, < r < 1, ip e [0, 27r):) 



-(t")=r" 

4r£>r(t) 



dTe iTE/n 



(1-r 2 ) 2 

r(t')=r' il>(t')=ip' 



x exp 



r 2 + r 2 -i/> 2 
2m- ^TTT + ft 



(1-r 2 ) 2 



— ?r / dp p sinh 7rp exp 
27r2 ^° ^ 



(1-r 2 ) 2 



32mr 2 



dt 



2m 



p 2 + 



x jW'-ii>')<p- 



l + r /z 



ip-f V 1 _ T i 



lM „/i 

2 

1 + r //S 
ip-f \ 1 — r //S 



+ ip + / 



(4.60) 



Here we have also stated the (equivalent) representation in terms of the Poincare disc via 
z = x\ + 1x2 = re 1 ^ = tanh ^(sin^> + icos</>) which is important is various discussions 
in the theory of quantum chaos [5-8, 61]. 
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4.3.5. General Pseudo-Ellipsoidal Coordinates. We consider the general elliptic 
coordinate system 



1 



Ur 



= , (i = 1,2) . 



(4.61) 



Pi - a pi - b Pi - c 
Explicitly in terms of the variables u on the pseudosphere A^ 2 ) they are given by 



«n = 



(Pi -C)(p 2 - C) 

(a — c)(b — c) 



Mi = 



(Pi ~&)(P2 -&) 

(a-6)(c-6) ' 



Mo = 



(Pi ~a)(P2 - a) 
(b — a)(c — a) 



(4.62) 

One now can make the identification p\ = —k' 2 nc 2 (p, k), p2 = k 2 nc 2 (z/, k'), a = —k', 
b = k 2 and c = 0. This yields 



Uq = nc(/U, k) nc(z/, fc') , /c 2 + /c' 2 = 1 , 
ui = dc(/U, fc) sc(z/, k') , 
it2 = sc(/U, k) dc(z/, fc') , 



(4.63) 



together with < k 2 < k 2 nc 2 a < —k' 2 < — k' 2 nc 2 a, say. The functions sc(p, k), 
dc(p, /c) and nc(p, k) denote Jacobi elliptic functions which correspond to the functions 
sn, dn and cn by taking the argument imaginary, where e.g. nc(/i, k) has the periods 
4\K, 2K + 2 \K', AK (K, K' the elliptic integrals) [1]. Note the relations nc 2 p + sc 2 p = 
1 and dc 2 p = 1 + k 2 sc 2 p. The metric tensor g a f, in these coordinates is given by 



{gab) = (k' 2 nc 2 p + k 2 nc 2 is)t 2 , 
and the momentum operators are 
h ( d ,2 sc/inc/idcp 



(4.64) 



P» 



— k 

i \dp k' 2 nc 2 p + k 2 nc 2 v 



Vv 



h ( d ^ 2 sc z/ nc z/ dc z/ 



i\du k' 2 nc 2 p + k 2 nc 2 v 



The Hamiltonian has the form 



(4.65) 



ft 



2m 



^A( 2 ) = 



d' 2 



<P_\ 

2m k' 2 nc 2 p + k 2 nc 2 u\dp 2 dv 2 ) 



■■(jpI+pI)- 



^ m a/ k' 2 nc 2 p + k 2 nc 2 v M " a/ /V 2 nc 2 p + k 2 nc 2 z/ (4.66) 
The path integral can therefore be formulated yielding 



ft ./o 



A.(t")=/i" 

dEe iET / h / P/i(t) 



i/(t")="" 



Vv(t)(k' 2 nc 2 p + k 2 nc 2 z/) 



2 2 



x exp 



At(t')=M' i/(t')=u' 
t" 

^/ (fc /2 nc 2 p + fc 2 nc 2 z/)(/i 2 + z> 2 )rft 
2n y t / 
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M(s")=M" 



v{a")=u" 



= ds" J Vfi(s) J Vu(s) 



M(0)=M' 



x exp <^ — 



v(Q)=v' 



ds 



dp 



^/o W(p* + \)/2m-E' 

^2 Q-i /2 -iV 2m^-i ( cosh d ^ (q". q')) 



(4.67) 
(4.68) 



where d A ( 2 ) (q", q') must be expressed in ellipsoidal coordinates. The eigen-functions 
^p,k(p>, ^) := A PjK (h)B PjK (v), are generalized hyperbolic Lame polynomials in ^ and v 

and are the normalized solution of the Hamiltonian — £— A A (2>: 

2m A ^ ; 



d 2 „4 P;K (/x) 
d^ 2 

d 2 ff p ,» 
dv 2 



+ 



+ 



| + (p 2 + i)^ 2 ncV 



+ (P 2 + i)k 2 nc 2 v 



Ai, K (fj) = , 
B hK {y) = . 



(4.69) 



Depending on the domain, the values p, and v can take on, one can discriminate six 
different coordinate system which are encoded in (4.62, 4.63). For this purpose one 
writes the kinetic energy-term in terms oi p\, 2 [69, 90] 



ill - ill- ill = \(Pi ~ Pi) 



p\ 



Pi 



P 3 (Pl) Ps(P2) 



(4.70) 



where the polynomial P 3 (p) is defined by Ps(p) = (p — a)(p — b)(p — c). Therefore we 
obtain the alternative path integral formulation 



Pi(t")=p" 



Pi(t')=p' 1 



P2(t")=P2 
P2(t')= P ' 2 



(-i)(pi - P2) 

4y/P 3 (p 1 )P 3 (p 2 ) 



x exp 



ft Jt> 

Pi{s")=p>> 



m 



(Pi ~ P2) 



P\ 



Pi 



Ps(Pl) Ps(P2) 

P2( S ")=P2 



-AV PF 



dt 



J ds" J v /-P 3 (p l )Vp l (s) J ^P, 



, 3(p 2 )Vp 2 (s) 



Pi(0)=p' 1 



P2(0)= P ' 2 



xexp ^nJo 



f (pMpI - pMpI) + - ^ £ m, (r? + 2r 

i=i 



2 D 



dp 



^Jo h 2 (p 2 + l)/2m-E 



a;j p [)b;jp' 2 )a p M)b p Ap2) 



m 



7rn 2^-l/2-i^/2mE/h 2 



-(coshrf A( 2)(q // ,q / )) , 

4 
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ds 

(4.71) 
(4.72) 



with the functions A PjR and B P)K rewritten in p\2 and cosh<i A (2) expressed in pi^. 
Furthermore denote M x = l/P 3 (p 2 ),M 2 = P 3 (pi), S = (pi - p 2 )/4P 3 (pi)^3(P2) and 
Ti = P^(pi)/2P 3 (pi) (i = 1, 2). The six different coordinate systems are now defined by 
by the roots of the polynomial Pz{p) (simple real, complex, double or triple real) with 
the requirement of the positive definiteness of the two quantities 

^—^ > , -£f^ > . (4.73) 

They are given by [69, 90] 
4.3.5.1. Elliptic Coordinates. 



c<b<p2<a<pi 



u{ ui u z 



+ — ^- = 0, (< = 1,2) 



Pi - b pi - a Pi - c 
4.3.5.2. Hyperbolic Coordinates. 

p 2 <c<b<a<pi 



Pi -a [(Pi - 7) 2 + 5 2 ] 

4.3.5.4. Elliptic-Parabolic Coordinates. 



Pi - a pi-b {pi- b) 

4.3.5.5. Hyperbolic-Parabolic Coordinates. 

p 2 <c = b<a<pi 



4.3.5.6. Semi-Circular-Parabolic Coordinates. 



p 2 <c = b = a<pi , 

U -U! \ 2 2 2 

h« 2 =U Q -U 1 , (l= 1,2) 

Pi - a J 



(4.74) 



« 3 + ^L_^L =0i (j = 1>2) . (4-75) 



Pi - c pi - a pi 

4.3.5.3. Semi-Hyperbolic Coordinates. 

P2 < a < pi , 6 = 7 + i<5, c = 7 - i <5 
W 2 2 < 5« 1 « + (Pi-7)(^-«j) =0 (z = i2) (4 - 76) 



c = 6<p 2 <a<pi , 

u| , ul-u\ {u - Ul ) 2 (4.77) 

= , (z = 1,2) . 



ul . ^-^ 2 (txo-m) 2 n 1 9 . (4-78) 

H r 7 rry = » (i = 1, 2) . 

i - a Pi - o (pi - 6) 2 



(4.79) 



This concludes the discussion of the separable coordinate systems on two-dimensional 
spaces of constant curvature. 
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5. Separation of Variables in Three Dimensions 

In this Section we give the path integral representations in the three-dimensional spaces 
of constant curvature E^ s \ and A^ 3 ). The notation and presentation will be as in 
te previous Section, with the only difference that the stating of the invariant distances 
will be given separately in each case. 

Some new path integral evaluations will be presented, however, our emphasize again 
is more concerned with the connection of the various path integral identities with each 
other. The wave-function expansions then represent again the various possible group 
matrix element representations yielding identities connecting the coordinate realizations. 
We also note the general form a potential must have to be separable in a coordinate 
system. The two-dimensional analogue then follows from omitting the ^-coordinate. 



5.1. The Flat Space 



5.1.1. General Form of the Propagator and the Green Function. Similarly as in the 
two-dimensional case, the propagator and the Green function on are well-known, 
and are given by 



K S(3) (dWq",q');) = 

g e(3) (d EW (^y,E) 



m 



3/2 



2n\hT 
m 



- ^4(3) (q",q') , (5.i) 



2\hT 



q",q') P V & J 



4nfi 2 d E(3) (q",q') \ h ) ^ 

q denotes any of the eleven coordinate systems which allows a separation of variables 
of the Laplacian in three dimensions. Due to the rather many coordinate systems, the 
explicit expression for d E (3) (q", q') will be given together with the definition of the 
coordinates. In the nomenclature we follow [69, 87]. 

5.1.2. Cartesian Coordinates. Again, we start with the simplest case, cartesian 
coordinates (x, y, z) = x e IR 3 . Then (g a b) = H 3 , 



4 (3) (q",q') = |x"-x'| 2 , 
and p = — i W. This gives for the Hamiltonian 



h 2 



2m^ E<3) 2m 
and for the path integral we have [31, 32] 

x(t")=x" 



2m' 



(5.3) 



(5.4) 



/ 



im 



Dx(t)exp( ^- / ± 2 dt\ = 



m 



x(t')=x' 



2n\TiT 
dp 



3/2 



exp 



im 



IR 3 (2tt) 



exp 



\2HT 



11 i\2 
x — X 



2m 



p 2 + i p • (x" - x') 



(5.5) 



(5.6) 



Note the formulation via a path integration over the Euclidean group [11]. The most 
general potential separable in cartesian coordinates has the form 



V = u(x) + v(y) + w(z) 



(5.7) 
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5.1.3. Circular Cylinder Coordinates. Next we consider circular cylinder coordinates 
which are very similar to the two-dimensional polar coordinates 



x = r cos <f> , r > , 

y = r sin , < <fi < 2% , 

z = z , z G IR . 



(5.8) 



Here we have for d E (3) (q", q') 



4(3) (q", q) = \z" -z'\ 2 + r' z + r" z - 2r'r" cos(0" - cf> f ) 



(5.9) 



The metric reads (g a b) — diag(l, r 2 , 1), and the momentum operators are given by (4.10) 
and p z = — i hd z . This gives for the Hamiltonian 



2m~ iE<3> 2m\dr 2 r Or r 2 d<fi 2 ' dz 2 J 2m\^ r ' r 2 
We therefore obtain the path integral identity [4, 10, 55, 92, 104] 

r(t")=r" 4>{t")=4>" z(t")=z" 



H 2 ld_ : 1 d 2 d 2 

+ -7T- + ~ n TTTn + 



8mr 2 
(5.10) 



y r^r(t) y v^t) y 

4>(t')=<f>' z(t')=z' 



r(t')=r' 



x exp 



(ft 



m \ ' 
2^iWf) 



3/2 



exp 



lm / 



2ftT 



(V'-/| 2 + r' 2 + r" 2 ) ^e^^Hf™ 



mr r 



hT 



(5.11) 



IR. 



^ 2 i 



p z (z"-z')-ihp 2 z T/2m 



E 

Ze^ 



A separable potential has te form 



?K<t>" -4>') 
2^ 



1 



pdpJi{pr')J l (pr")e- ihp2T/2rn . 

(5.12) 



F = u(r) + -^v{4>) + w(z) . 



(5.13) 



5.1.4. Elliptic Cylinder Coordinates. We consider the coordinate system 
x = d cosh \i cos v , ^ > , 

y = <i sinh \i sin z/ , — n < is < n > (5-14) 

z = z , z G IR . J 

(alternatively G IR, < z/ < 7r [84]). Here the invariant distance is 

d E ( 3 ) (q", q') = 1 2" — z'| 2 + ^(cosh^' cosily" cos z/ cos v" + sinhju' sinhju" sin is' sin v") . 

(5.15) 
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The metric is (g a b) = diag[(i 2 (sinh 2 /i + sin 2 z/), d 2 (sinh 2 \i + sin 2 f),l], and we obtain for 
the momentum operators (4.16), and p z = — ift<9 2 . Consequently for the Hamiltonian 



2m 



2m 
1 
2m 



/ d 2 



d 2 



d 2 



d 2 (sinh 2 ^ + sin 2 z/) \ d[i 2 du 2 ) dz 2 
^(sinh 2 ^ + sin 2 is)~ 1/2 (pl + pi) (sinh 2 fi + sin 2 z/)" 1/2 + p 2 



d' 2 



(5.16) 



The path integral construction is straightforward, however, again no explicit path 
integration is possible for the coordinates fx and v. The corresponding wave-functions 
yield again as eigen-functions Mathieu functions. Because we know the kernel in 
in terms of the invariant distance d E (3) we can state the following path integral identity 



M(i")=M" 



v{t")=v 
jl{t')=jl' v(t') = v> 



■,")=v" z(t")=z" 

J Pz/(t)rf 2 (sinh 2 ^ + sin 2 z/) J Vz(t) 



z(t')=z' 



x exp 



i m I 

2hJ t , 



d (sinh \i + sin z/)(/i + v ) + i 



.•2 



(ft 



dpz i 



p z ( 2 "-2')-i?ip^T/2m 



27T 



2tt ^ 



4 



p<ip e" 



ihp 2 T/2m 







m 



x meKV, ^-) me„fa", ^-) Me« *(£', Me«(£ 

3/2 



' 4 



27ri/iT 



exp 



i m 
2hT 



4o) (q",q') 



(5.17) 
(5.18) 



and d_E(3) must be taken in elliptic cylinder coordinates, and we have used the same 
notation as for the two-dimensional elliptic coordinates. A potential separable in these 
coordinates reads 

. (cosh + , (cos ,) + ^ (5i9) 
sinh /j, + sin v 



V = 



5.1.5. Parabolic Cylinder Coordinates. The last example for cylinder coordinates in 
three dimensions are the parabolic cylinder coordinates 



y 



\{r, 2 -e) , £GR,»7>0 , 



z = z , z e H , 

which is the obvious generalization of the two-dimensional case. Therefore 



(5.20) 



d|(»>(q'W) 



I II 1 1 2 
\Z — Z \ 



+ 



1 



-//2n2 



//2 



J„H tl til 



(5.21) 



We have (g^) = diag(£ 2 + ?? 2 , £ 2 + ?? 2 , 1), hence for the momentum operators (4.23) and 
p z = — ihd z . This gives for the Hamiltonian 



2m 



A_E(3) = 



2m 



1 



a 2 a 2 
+ 



£2 _|_ ^2 ^5^2 ^2 

30 



+ 



dz 2 



2m 



(e+v 2 r 1/2 (Pi+P 2 v )(e+v 2 )- 1/2 +P 2 z 



(5.22) 



The arising path integral in these coordinates must be treated similarly as in the two- 
dimensional case, the only difference being the z-path integration which is trivial. We 
obtain 



£(*")=£'■ 



v(t")=v' 



z(t )—Z 



J mt) J v V ( t )(e+v 2 ) J vz{t) 



x exp 



i m r f 

2h J t , 



z{t')=z' 

(e+v 2 )d 2 +v 2 )+z 2 



dt 



IR 



d Pz e ip z (z"-z')-ihp 2 z T/2m 

2tt 



e,o ^ ^ 

in the notation of 4.1.5. A separable potential has the form 

u(0 + v(v) 



V = 



e+v 2 



+ v(z) . 



(5.23) 



(5.24) 



5.1.6. Sphero-Conical Coordinates. We consider the coordinate system 

x = r sn(/x, k) dn(i/, k') , r>0 , 
y = r cn((U, k) cn(z/, k') , k 2 + k' 2 = 1 , 
^ = r dn(jU, k) sn(z/, fe') , 

in the notation of 4.2.3. The invariant distance is given by 

4o) (q", q') = r' 2 + r" 2 - 2r'r" ( sn fj," sn //' dn z/' dn !/' 

+ cn fx" cn cn i/" cn v' + dn dn \s! sn i/' sn v'^j 

The metric tensor g a 6 in these coordinates has the form 

(gab) = diag [l, r 2 (k 2 cn 2 \i + k' 2 cn 2 u),r 2 (k 2 cn 2 ^ + k' 2 cn 2 z/)] , 



(5.25) 



(5.26) 



(5.27) 



and the momentum operators are p r = — i h(d/dr + 1/r), and p„ as in (4.40). The 
path integral can now be formulated yielding 



r(t")=r'' 



/*(*")=**" 



u(t")=u' 



J r 2 Vr(t) J Vp.it) J Vv(t)(k 2 cn 2 /j + k' 2 cn 2 v) 



r(t')=r' 



M (t')=M' 



v(t')=v' 
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x exp 



t" 

^ / (r 2 + r 2 (k 2 cn 2 ^ + k' 2 cn 2 u)(ft 2 + u 2 )) 
2a J t i V / 



dt 



r(t")=r" 



1 = K 



X 



Vr{t) exp 



r(t')=r' 



ft /„ V 2 



2m 



m 



iftTVrV" 



exp 



m / /2 . //2s 

(r + r ) 



2iftT 



x E E A l %(Ai / )5 l %(^)^,«M5«,«MJ« + i/2 ( - 



K 1=0 



mr'r" 



TiT 



(5.28) 



« fc£ Jo ^ r r (5.29) 



2tt 



m \ 3/2 

i^rj exp 



m 



2\hT E(3) 



(5.30) 



with d S (3) (q", q') in sphero-conical coordinates. Here a separable potential must have 
the form 



r 2 & 2 cn 2 a + k' cn 2 /? 



5.1.7. Spherical Coordinates. We consider the spherical coordinates 

x = r sin 9 cos , r > , 

j/ = r sin (9 sin ^ , < 9 < n , > (5.32) 

2 = r cos 9 , < < 2/T . 

These are the usual three-dimensional polar coordinates. Here 

4o) (q", q') = r' 2 + r" 2 - 2r'r" ( cos 0' cos 0" - sin 9' sin 0" cos(0" - 0')) • (5-33) 

The metric tensor is (<7 a b) = diag(l, r 2 , r 2 sin 2 6*), and the momentum operators have 
the form p r = — i %{d/dr + 1/r) together with (4.33). For the Hamiltonian we obtain 



2m 



A_E(3) = 



2m 



a 2 2 a i / a 2 3\ i 

+ 7T~ + ~o 7^7 + COt0 — + 



^^y»2 y» ^^f* y»2 I ^^^2 

,2 



00 



5 



2 1 



2 sin 2 d</> 2 



1/2 1 2 1 

2m \ r A 



r 2 sin 9 



8mr 2 



1 + 



1 



sin 2 



(5.34) 



The corresponding path integral is well-known [22, 37, 92] and we have the identity 

r {t")=r" e(t")=e" 4>{t")=4>" 

J r 2 Vr{t) J sin0£>0(£) J V<f>(t) 



r(t')=r' 



0{t')=0> 



4>{t')=<l>' 



32 



ft" V 
~h j t IT + r2 ^ 2 + r2 sin2 ^ 2 ) + 



1 + 



1 



8mr 2 V sin 2 6> 



dt 



m 



An i hTVrY 7 

00 

x 

z=o 

00 i 



exp 



171 / /2 . //2s 

(r + r ) 



2iftT 



X)(2/ + m ( cos (q", q / ))^+i/2 



^ » 1 /-00 

= E E YT*(e , 1 <l/)Yr(6" i r)-n= / pdpJj+iCprOJ.+ iCpr'Oe 

A potential which is separable in spherical coordinates reads 

V = u{r) + + -^-^-^(0) • 

r z r z sin 



5.1.8. Parabolic Coordinates. We consider the coordinate system 



x = £?ycos , 



y = C?7sin0 , f, 77 > , 

0< </>< 2tt . 



This gives for the metric tensor (g a b) = diag(£ 2 + ?? 2 , £ 2 + rj 2 , £ 2 ?y 2 ), and 



4 (3) (q", q ') = 7 (r+n 2 +(/+e' 



2(r/ 2 - r 2 )(V /2 " r' 2 ) - cos(0" - 0') 



For the momentum operators we get 



h ( d 



= 7 



+ 



+ 



1 



ft / d 



together with = — i ft<9<£. For the Hamiltonian we obtain 



+ 



V 



2m 



2m 
1 
2m 



d 2 d d 2 d\ Id 2 



£2 _|_ ^2 ^5^2 ^2 ^ y £2^2 £02 



/ 2 , 2\ 1 1 2 

(^+P^)^== + 7zl2^ 



8m^ 2 r] 2 



We obtain the path integral identity [15, 50] 

«(*")=€" »?(*")=»/" 



0(t')=4 
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X 



+ 



€(*")=€" 

e -il{(j>" -$') 



= (cwr 1/2 £ 



8m£ 2 ?y 2 
»»(*")=»»' 



2tt 



XeXP H./, 



£(*')=€' »»(*')=*»' 

2/2 1 l 



f« 2 + ^ + -> 2 )-^W 



eft 



hp 2 T/2m 



2t V Jo P 47r2eeWT 4 (l + 



= E 

x M_ iC/2p ,| Z | /2 (-ipr 2 )M iC/2p ^| /2 ^^ , 

(5.42) 

where the last line was obtained in a similar way as the corresponding result in 4.1.5, 
and use has been made of (A. 2. 2). In parabolic coordinates a potential is separable if 
it has the form 



y = HO +V(V) W{4>) 



^x 2 + y 2 + z 2 ' x 2 + y 2 



(5.43) 



5.1.9. Prolate Spheroidal Coordinates. We consider the coordinate system 
x = d sinh fi sin v cos <fi , fx > , 

y = d sinh /i sin v sin </> , < z/ < 7r , > (5.44) 

2 = <i cosh fi cos z/ , < <j) < 2n . ) 

*~y 2 2 2 2 2 2 

This yields (g a b) = d diag(sinh fx + sin z/, sinh + sin z/, sinh jusin z/), and for the 
momentum operators we obtain 



hid sinh u cosh a 1 \ Hid sin v cos z/ 1 

^a* = ~ 7T + ~r Tl — ; — r~2 — H tt coth A* h ^ = T I ~ + ~r — — : — — — ^ cot ^ h 



i \ sinh 2 u + sin 2 z/ 2 



i \ dv sinh 2 u + sin 2 v 2 



and p^, = — i hd^. The invariant distance is given by 



(5.45) 



d 2 E(3) (q", q') = d 2 cosh fi' cosh ju" cos v' cos z/'+sinh // sinh fi" sin z/ sin z/' cos(<//'— </>') 



The Hamiltonian has the form 



(5.46) 



^ A 

h 2 
2md 2 
1 



a 2 cQth d d 2 d \ 1 d 2 

sinh 2 + sin 2 v \ dfi 2 C ° ^ 9// <9z/ 2 °° " dv) sinh 2 ^ sin 2 z/ d 2 4> 



2md 2 



1 



V sinh /X + sin z/ 



+ 



1 



a/ sinh 2 + sin 2 v sinh 2 |U + sin 2 v 



fii 



h 2 



8md 2 sinh \x sin f 



(5.47) 
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The path integral construction is straightforward, however no explicit path integration 
is possible. Actually, an expansion into the corresponding wave-functions in the 

coordinates fx and v yields the spheroidal functions ps™(cosz/, 7) and S™ (cosh /x, 7) 
(7 2 = 2mEd 2 /h 2 = p 2 d 2 ), respectively, as the eigen-function of the Hamiltonian, a 
specific class of higher transcendental functions [84] , similar to the Mathieu functions. 
However, because we know on the one side the eigen-functions of the Hamiltonian in 
terms of these functions [84], and on the other the kernel in in terms of the invariant 
distance d E (2) we can state the following path integral identity (note the implemented 
time-transformation) 



a»(*")=m" 



4>{t")=4>' 



J V/j(t) J Pz/(t) d 3 (sinh 2 ,u + sin 2 z/) sinh ,u sin z/ J V<f>(t) 



XeXP ^, f 



m 



4>{t')=4>' 

2 . -2 ..12 ' 



'-^■d 2 (Jsinh 2 [i + sin 2 v) (/i 2 + z> 2 ) + sinh 2 [i sin 2 z/</> 2 j 



+ 



8md 2 sinh fi sin v 



dt 



= (sinh fj,' sinh /j," sin v' sin u") x l 2 ^ 



2^ 



x 



^ -iET/fi [°° ds" 



IR 



2tt/i 



xexp^- 



i/(s")=i'" 

p/i(s) y x>^(s) 

H(0)=n' v(0)=v' 



ft 2 / / 2 - i 

m ' -.2 , ,-.2\ , rJ 2/„;„u2 .. , „;„2..\ 1(1 4 



— (A 2 + z> 2 ) + £d 2 (sinh 2 /I + sin 2 i/) -— . 
Z Zm \sinh sin v 



+ 



I 2 - i 

1 4 




2tt 



neA 



p 2 dp 



2n+l(n-iy._ ihv * 



(n + l)\ 



hp 2 T/2m 



x ps^ ( cos v' , pd) ps^ ( cos v" , pd) S 1 ^ 1 } * ( cosh // , pd) 1 } * ( cosh //' , pd) 

3/2 



m 



2yri/iT 



exp 



lm 

2nr 



4o) (q",q') 



(5.49) 
(5.50) 



and d S (3)(q",q') must be taken in prolate spheroidal coordinates. The functions 

ps^(cosz/, 7) and Sh 1 (cosh /x, 7) are mutually determined by the separation parameter 
A = A^(7) giving an infinite countable set {A n }, (n e A). The functions ps^(cosz/, 7) 
yield in the limit 7^0 the associated Legendre-polynomials P^(cosz/), and the func- 
tions St (cosh/x, 7) the spherical Bessel functions, i.e. S l n l \z/^,^) oc y / iv/2zJ l+ i (z) 

(7 — > 0), therefore obeying the correct boundary-conditions. A separable potential must 
have the form 



u(cosh fj,) + v(cos v) 



sinh 2 \i + sin 2 z/ sinh 2 \i sin 2 1/ 



(5.51) 



Note that by the replacement z — > d(cosh / ucosz/ + 1) we obtain the spheroidal 
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coordinate system in which the Coulomb-problem in IR 3 is separable [78], e.g. c.f. [50] 
and references therein. 



5.1.10. Oblate Spheroidal Coordinates. We consider the coordinate system 



x = d cosh £ sin v cos < 
y = d cosh £ sin v sin < 
z = d sinh £ cos v , 



£>0 , 

< V < TV , 

< 6 < 2tt 



(5.52) 



(alternatively |i G K, < y < ix/2 [84]). This yields (<7ob) = d 2 diag(cosh 2 £ — 

2 2 2 2 2 \ 

sin u, cosh £ — sin v, cosh £ sin v) , and for the momentum operators we obtain 



hid sinh £ cosh £ 1 



+ 



i V d(i cosh 2 £ - sin 2 v 2 



+ - tanh £ , p v 



h / d 



+ 



sin v cos z/ 



i I dv cosh 2 £ - sin 2 v 2 



1 

+ - cot V 



and = — i ftc^. The invariant distance is given by 



(5.53) 



d 2 E(3) (q", q') = d 2 sinh £' sinh £" cos z/ cos i/' + cosh £' cosh £" sin z/ sin i/' cos(0" -</>')] • 

(5.54) 

The Hamiltonian has the form 



^ A 

h 2 



2md 2 
1 



1 



2md 2 



cosh 2 £ - sin 2 v \d£ 2 

/ ,2 C ■ 2 = (Pt+Pt 
V cosh £ — sin z/ 



a 2 , , a a 2 d , 

+ tanh£ — + + cot z/-^— ) + 



<9£ <9z/ 2 
1 



1 



d 2 



+ 



dv ) cosh 2 £ sin 2 u d 2 (j) 
1 



A/cosh 2 £ - sin 2 u cosh £ - sin v 



-P<t> 



8md 2 cosh £ sin f 



(5.55) 



The path integral construction is straightforward, however again no explicit path 
integration is possible. An expansion into the corresponding wave-functions in the 
coordinates £ and v yields eigen-function of the Hamiltonian in terms of spheroidal 

functions ps™(cos u, ij) = psi™(cos u, 7) and S , ^T ( ' 1 ' ) (— icosh£, i7) = Si™ ^(cosh^, 7) 
(7 2 = 2mEd 2 /h 2 = pd), respectively. Because we know on the one side the eigen- 
functions of the Hamiltonian in terms of these functions [84], and on the other the 
kernel in in terms of the invariant distance d E (3) we can state the following path 
integral identity 



«*")=€" 



<t>(t")=<t>" 



x exp 



v(t')=u 

t" 



h 



m 



2 ,-„;„2 ,.12 



d 2 (coslr i - sin 2 v) (£ 2 + z> 2 ) + coslr £ sin 2 v 



+ 



8md 2 cosh £ sin v 



dt 
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= (cosh £' cosh £" sin v' sin z/') 

lezz 



e ii(^"-^') 
2^ 



i/(«")=i/" 



x exp 



27r^ 
i 

n 



J ds" J V^s) I Vv{s) 



rn 



u(Q)=u' 

2 „;__2 



(r + v ) + Ed 2 (cosh 2 £ - sin 2 i/) 



ft 2 / / 2 - i / 2 - \ 



- 



2m \ cosh £ sin z/ 



e il(*"-^') 2 

^ 27T ^ 7T 

n6A 



' ^ 2n + l (n-Q! p _ s ^ T/2m 



(n + Z) 



x p< (cos z/, pd) psi^ (cos i/", pd) Si^ (1) * (cosh pd) Si^ (1) * (cosh f ", pd) 

3/2 



m 



2yri/T,T 



exp 



lm 
2hT 



4(.)(q".«f) 



(is 
(5.56) 

(5.57) 
(5.58) 



and <i S (3) (q", q') must be taken in oblate spheroidal coordinates. Due to the con- 
struction of the spectral expansion the functions psi^(cos v, 7) and Si^ 1 ^ (cosh /i, 7) are 
obeying the correct boundary-conditions. Here a separable potential must have the 
form 

W{(P) (5.59) 



'u(sinh^) +v(cosv) 



2 2 2 2 

cosh /i — sin z/ cosh ^ sin z/ 



5.1.11. Ellipsoidal Coordinates. The last two coordinate system are the most compli- 
cated ones and are similar in some of their features, c.f. Section 3.1. First we consider 
the coordinate system 



' " w a 2 (a 2 -b 2 ) ' y y 6 2 (6 2 -a 2 ) 

* = ^ , £ 2 > a 2 > £ 2 > 6 2 > £ 2 > c 2 = . 

The metric tensor is given by 

(a -am-® (a - am - e 3 ) (a -am-® 



(5.60) 



( ^ d "H K(S) v ■ m ) • (6 ' 61) 

with P 2 (£ 2 ) = (£ 2 - a 2 )(£ 2 - b 2 ). With the identification £1 = adn(A, fe)/cn(A, fc), 
£2 = adn((U, fc') and £3 = 6sn(z/, fc) (6 = ka, V« 2 — & 2 = fc'a = <i) this can be rewritten 
into 



sn(A, k) sn((U, k') dn(z/, fc) j 011 ^ cn(zv, fc) 

cn(A, cn(A, k) 

dn(A, fc) dn(//, fc') sn(z/, fc) 
cn(A, fc) 
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(5.62) 



Due to the very complicated structure it is of now use to write down all subsequent 
necessary quantities for the path integral. Instead we exploit the results of Section 3.2, 
in particular the separation formula (3.21). We identify (g a b) = diag(/if , h 2 ,, h 2 ), 
furthermore 



r, : = 



(i = 1,2,3) 



' 

1 l/(a-a 2 ) l/(a-b 2 )(a 2 -b 2 ) 

1 1/(^1 -a 2 ) l/(^ 2 2 -6 2 )(a 2 -6 2 ) 

1 l/fe 2 -a 2 ) l/(il-b 2 )(a 2 -9) 



(5.63) 
(5.64) 



Mi = 
M 2 = 
M 3 = 



2 - &h (e 2 2 - a 2 ) (el - (e 2 - v) (ei - « 2 ) 



& 2 V(ei 2 -& 2 )(el-a 2 ) (e 2 -a 2 )(ef-& 2 ) 



1 



1 



1 



& 2 V(e 2 -a 2 )(e 2 2 -& 2 ) (ei 2 -& 2 )(e 2 2 -« 2 ) 



(5.65) 



and obtain the following path integral identity 
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eft 
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x exp 
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j 8m 



^M,(r 2 + 2r:) 
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./o 
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(5.66) 
(5.67) 

(5.68) 

and d E (3) (q", q') must be expressed in ellipsoidal coordinates via rf S (3) = |x({e"} — 
x ({e'})l- The functions * /s ,a, p (Ci,6»C3) = ^Ap^O^Ap^C^A^fe) are solutions 
of a three-fold Lame equation with separation parameters X,K,p, respectively [87, 
p. 1305]. A potential separable in ellipsoidal coordinates must have the form 



3 /2 

4(3) (q",q') 



2yriftT 



exp 



2iftT 



V = 



(e| - e 3 2 Mei) + (e| - e 3 2 Me 2 ) + (e 2 - ejMegj 
(ei 2 -e 2 2 )(ei 2 -el)(e 2 2 -el) 



(5.69) 
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5.1.12. Paraboloidal Coordinates. As the last coordinate system in we now consider 



x = 



(£i 2 -a 2 )(£ 2 2 -a 2 )(£f-a 2 ) 



(a 2 - b 2 ) 



(b 2 - a 2 ) 



z = ^{ei+e2+e 3 -a 2 -b 2 ) , ei>a 2 >e 2 >b 2 >es>c 2 







The metric tensor is given by 

(„ \ A\~* f c 2 ^1 ~~ £ 2 )(£i ~~ C3) ^2 (^2 €3) £2 (£3 — ~~ £2) 

\9ab) = diag £i ——27 ,f 2 ^772n »?3 — 



(5.70) 



mi) 



Pirn 



Pirn 



(5.71) 

Similarly as before we can make the identification £1 = adn(A, fc)/cn(A, k), £2 = 

adn(z/, k') and £3 = a^l - k' 2 /cn 2 (//, fc) (6 = fca, Va 2 - 6 2 = fc'a = v 7 ^) and (5.70) 
can be rewritten into 



x = d 



sn(A, fc) sn(z/, fc') 
cn(A, fc) cn( / u, fc) 



y = d 



sn(/x, fc) cn(^, /c') 
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sn 2 (A,fc) sn 2 (^,A;) dn 2 (i/,fc') 



cn 2 (A, fc) cn 2 (/x, fc) 



fc' 2 



(5.72) 



We proceed similarly as before and use (3.21) for the path integral formulation. We 
identify (g a b) = diag(/i 2 , /?|, /i|), furthermore 
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el e 3 2 /(el-a 2 ) ^/(el-W 



6 2 ) 
6 2 ) 
6 2 ) 



(5.73) 
(5.74) 
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e?e 3 2 / 1 1 
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Ms = a^ 2 " V(ei 2 -« 2 )(e 2 2 -^ 2 ) (e! 2 -6 2 )(ei-a 2 ) 

and obtain the following path integral identity 
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= (s's'T 1/4 I 
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(5.76) 
(5.77) 

(5.78) 



and g?e(3) (q", q') must be expressed in paraboloidal coordinates via <i E (3) = |x({£ ;/ } — 
x({f'})|. The functions # k ,a, p (£i, 6, 6) = -4 k ,a,p(6)5 k ,a,p(6)C k ,a, p (6) are solutions 
of a three-fold Lame equation with separation parameters \n,p [87, p. 1305]. A 
potential separable in ellipsoidal coordinates has the form of (5.69). 



5.2. The Sphere 



5.2.1. General Form of the Propagator and the Green Function. We consider the sphere 
as the other space of constant positive curvature. There exist six coordinate systems 
which admit separation of variables on which will be discussed in the following. 
The propagator and the Green function on are best calculated in terms of polar 
coordinates. One obtains [12, 59, 99] 



(5.79) 



rs™(,i, ( n " n '\-jA m sin [( 7r -^ ) (q // ,q / ))(a+|) J 

G (Vs(3)(q ,q),-fc ) = — — 2— r~ - f — , ui ■ — ; Tli — a -- 

V / 2-k% sm[7r(a + f )J sin^ S (3)(q", q') 



where a = — \ + \J 2mE/h 2 + 1, and 6s(z\t) is a Jacobi theta- function [41, p. 931] 

oo 

9 3 (u,q) = 6 3 {u\t) = 1 + 2j2<l n2 cos(2™) , (q = e i7rr ) . (5.81) 



n=l 



5.2.2. Spherical Cylinder Coordinates. We first consider the coordinate system 

so = cos 9 cos 0i , < < 7r/2 , 
si = cos flsin^! , 

s 2 = sin cos 02 , < 0i 5 2 < 27T , 
S3 = sin 6* sin </> 2 . 
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The metric reads (g a b) = diag(l, cos 2 0, sin 2 9), and the invariant distance is given by 

cosV> S (3)(q",q') = cos 9' cos 0" cos« -</>[)+ sin 9' sin 9" cos^ - </>' 2 ) . (5.83) 
The momentum operators are 

ft / d 
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i V <30 



and p^ 12 = — ihdfa z. Therefore for the Hamiltonian 
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We obtain for the path integral 
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(5.86) 



5.2.3. Sphero- Conical Coordinates. We consider the coordinate system 



50 = cos , 

51 = sin0 sn((U, fc) dn(v, k') 

52 = sin cn(//, fc) cn(z/, A;') 

53 = sin0 dn(/i, fc) sn(z/, fc') 



tt/2 < < tt/2 , 



A; 2 + A; 



/2 



(5.87) 
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in the notation of 4.2.3. The metric has the form (g a b) = diag[l, sin 2 9(k 2 cn 2 /j + 
k' cn 2 z/), sin 2 9(k 2 cn 2 (i + k' 2 cn 2 u)] and 

cos i/j S (3) (q", q') = cos 9' cos 9" — sin 9' sin 9" 
x sn fx" sn // dn z/' dn v' + cn ju" cn \s! cn z/' cn v' + dn //" dn // sn z/' sn z/ . (5.88) 

For the momentum operators we have 



Pe 



h ( d 



+ cot0 ] 



i \d9 j 

together with and p v as in (4.40). We have for the Hamiltonian 

1 



(5.89) 



/I 2 n 2 

~ 2m" As(3) " ~2m" 



d 2 n n d 
-2cot0 — + 



1 

2m 



P 2 e + 



89 2 
1 



09 sin 2 9(k 2 cn 2 /j + k' 2 cn 2 z/) V <V <^ 2 



a 2 a 2 

+ 



pI+pI 



1 



aA 2 cn 2 + /c' 2 cn 2 z/ sin a/ /c 2 cn 2 \x + k' 2 cn 2 z/ / 2m 



(5.90) 



The path integral can therefore be written down yielding 



e(t")=e' 



/i(t")=M" 



i/(t")=f" 



y sin 2 0£>0(t) y £>//(£) J Vv(t)(k 2 cn 2 ij + k' 2 cn 2 v) 

6" /j(t')=M' i/(t')=^' 

"~— (§ 2 + sin 2 0(/c 2 cn 2 [i + k' 2 cn 2 z/)(/i 2 + z> 2 )"j + ^— 



6(t>)=9> 



XeXP ^/. 



(sin (9' sin^)- 1/2 e" i?lT/2m £ £ m^W) 
6(t")=6" 



K 1 = 



X 



X>0(£) exp 



6{t')=6' 



h Ih 



h 2 (i + \) 2 - 1 



2m cos 2 



+ ^ + 2 ) 



K Z = 

oo 



E + *) (/- O! ) <^(sinn<v / 2 2 (sm0O exp 



inr 

~2m 



N{N + 2) 



(5.91) 



in the notation of 4.2.3, and we have set iV = n + /, n G 1N . 



5.2.4. Spherical Coordinates. We consider the usual three-dimensional polar coordinates 

50 = COS 01 , < 01,2 < 7T , 

51 = sin 0i sin 2 sin , < < 27r , 
s 2 = sin 0i sin 2 cos , 
s 3 = sin 0i cos 2 • 



(5.92) 
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Here we have (g a b) = diag(l, sin 2 0\, sin 2 6>i sin 2 9 2 ) and 

cos ^ S (3) (q", q') = cos 6[ cos 6" + sin 6[ sin 9'[ ( cos 6> 2 cos 6> 2 ' + sin 6> 2 sin 9 2 cos(0" - 0')) . 

(5.93) 

For the momentum operators we have (5.89) for po 1 , and for pg 2 and p<f> as in (4.33). 
For the Hamiltonian we obtain 



Ti h 

2^ As(3) ~ 2^ 



1 v d 2 8 1 / d 2 8 Id 2 

+ 2cot6>i— + — ^ — — + cot/9 2 — + 



1 



dOi L d9 x sin 2 6»i \d0 2 

.2 \ + 2 



89 2 ' sin 2 6» 2 ^ 



1 



+ 



1 



2m V 1 sin 2 6>i sin 2 6>isin 2 6> 2 / 8m V sin 2 9 X ' sin 2 9 X sin 2 6> 2 



(5.94) 



This gives the path integral formulation 

0i(t")=0'i e 2 (t")=^' 4>{t")=4>" 

J sm 2 9 1 V9 1 (t) J sin 9 2 V9 2 (t) J V</>(t) 

8i(t')=o[ e 2 {t')=e' 2 <f>(t')=<j>' 

t' : 



x exp 



h L/ 



in 



9\ + sin z 9 x 9l + sin 2 #i ^ ^ 



,2 a 62 



2 a 12 



h 2 / l l 

H 4 H =— H = = — 

8m V sin 2 9 sin 2 6>i sin 1 9 2 



^ £(' + Wi ( cos ^5(3) (q", q')) exp 



1=0 

oo 



ihT 
~2m 



1(1 + 2) 



E E *& 2 (*i , 6K 2 , <f>')*?Z,m 2 (*!', 91 <P") exp 

mi ,m 2 Z=0 



iftT 
~2m~ 



1(1 + 2) 



dt 



(5.95) 



(5.96) 



where the wave-functions are given by [28, p. 240] 

(5.97a) 



iV 5 (3) = 



2^32 — 1 — 27711 — 2m 2 



T(/ + mi + 2)r(mi + m 2 + 1) 



(/ + l)(mi + § (/ - mi)!(mi - m 2 )\ Y 2 ( mi + l)T 2 (m 2 + §) 



(5.97b) 



5.2.5. Twofold Confocal Ellipsoidal Coordinates 1). We consider the two-fold confocal 
ellipsoidal coordinates defined by 



4 + 4 

Pi-b 



+ 



+ 



'0 



Pi - a pi - c 



, (i = 1,2) , 



(5.98) 



43 



where c < P2 < b < p\ < a. The corresponding metric reads 

- (pi -b)(j>2-b)pl , 



pi 



pI 



ds 2 _1 _ , 

^ - 4 (pi ^vp 3 (pi) p 3 (p 2 ; 



(5.99) 



where P^(p) = (p—a)(p—b)(p—c). Using the notations of the sphero-conical coordinates 
we can identify 



s = sn( / u, k) dn(z/, k') , k 2 + k' 2 = 1 , 

51 = dn( / u, fc) sn(z^, k') cos , < <fi < 2n , 

52 = dn(//, /c) sn(z/, A;') sin , 

53 = cn(p,, fc) cn(z/, k') . 



(5.100) 



We omit details in the following. Obviously, these coordinates correspond to sphero- 
conical coordinates with an additional circular coordinate. This gives the path integral 
formulation 



H(t")=n" v{t") = u" <t>(t")=<t>" 

Vp(t) J Vv(t)(k 2 cn 2 p + k' 2 cn 2 v) V-k' 2 dn p sn p J T><f>(t) 



M (t')=M' 



x exp < — 



u{t')=u' 



<t>(t>)=4>' 



h 

e i %T/2m 



y [{k 2 cn 2 p + k' 2 cn 2 u) (p 2 + v 2 ) - k' 2 dn 2 p sn 2 za/> 2 ) - AV PF 
V's(3)(q",q / ) 



dt 



An 2 d cos ips( 3 ' 



■03 



TiT \ 
27rm / 



(5.101) 



and ips( 3 ) (q"> q') nas to be taken in these coordinates. Note that the (fi-path integration 
can be explicitly done, leading to a path integral in p, v with similar to the sphero-conical 
one with additional 1/ dn 2 p and 1/ sn 2 z/-terms, respectively. 



5.2.6. Twofold Confocal Ellipsoidal Coordinates 2). We consider the two-fold confocal 
ellipsoidal coordinates defined by 



+ 



s 2 + s 2 



+ 



, (i = 1,2) , 



Pi-b Pi - a Pi - c 
where c < pi < b < p\ < a. The corresponding metric reads 



(5.102) 



pi 



p\ 



ds2 _ 1 



+ (pi - a)(p2 - o)f>l 



(5.103) 



There is only little difference to the previous case. Proceeding similarly as before we 
have 



so = sn(p, k) dn(z/, k') cos <p , k 2 + k' 2 = 1 , 

51 = sn( / u, k) dn(z/, fc') sin0 , < 4> < 2tv , 

52 = dn(/z, k) sn(z/, fc') , 

53 = cn((U, fc) cn(z/, fc') . 



(5.104) 
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We omit details in the following. This gives the path integral formulation 

fi(t")=A»" u(t")=u" 4>{t")=4>" 



J Vp(t) J Vu(t)(k 2 cn 2 p + k' 2 cn 2 v) Vk 2 snp dn^ J D(j)(t) 



**(*')=/*' 



x exp H y t , 



v{t')=v> 



<t>(t')=<t>' 



y ((/c 2 cnV + k' 2 cn 2 v)(p 2 + z> 2 ) + A; 2 sn 2 ^ dn 2 z/0 2 ) - .A 5 



PF 



e i?iT/2m d . (q",q' 

-0 3 



47T 2 d COS ^5(3) 



27rm / 



(5.105) 



and V's(3)( t l"> q') nas to be taken in these coordinates. The difference of the last two 
coordinate system lies in the orientation with respect to an axis. 

5.2.7. Confocal Ellipsoidal Coordinates. Let us consider the coordinate system defined 
by 



+ 



+ 



0, (i = 1,2,3) 



Pi - d pi - c Pi - a Pi-b 
where d<ps<c<p2<b<pi<a (c.f. Section 3.1). Explicitly 



(5.106) 



(Pi ~ d)(p 2 - d)(p 3 - d) 
(a — d)(b — d) (c — d) 



Si = 



( Pi ~ c)(p 2 ~ c)(p 3 - c) 
(a — c)(b — c)(d — c) 



2 _ (pi - a)(p 2 - a)(ps - a) 2 _ (Pi ~ b)(p 2 - b)(p 3 - b) 
2 (d-a)(c-a)(b-a) ' 3 (d - b)(c - b)(a - b) ' 

The corresponding line element is given by 



(5.107) 



dt 2 



1 

4 



(Pl ~ P2){pl - Ps) .2 , (P2 ~ Ps)(P2 ~ Pi) . 2 . (pj ~ Pl)(P3 ~ gg) -2 

— Pi H ^77T^ Pi H ^77T^ P3 



P 4 (P1) rl p 4 (p 2 ) r/ A(p 3 ) 

(5.108) 

where -P 4 (p) = (p — a)(p — b)(p — c)(p — d) . We apply the separation formula (3.21). We 
identify (g a b) = diag(/if , /i 2 ., h 2 ), furthermore 



s -- 



(Pl - P2)(P1 - P3)(P2 - P3) 



Mi = 



4P 4 (P1)P4(P2)P 4 (P2) 
P2 - P3 



P 4 (P2)P 4 (P3) 



Mo 



P3 - Pl 



P 4 (pi)P 4 (p 3 ) 



Ma 



Pl -P2 



P 4 (P1) J P 4 (P2) 



(5.109) 
(5.110) 

(5.111) 



and obtain the following path integral identity 

Pi(*")=Pi 



f[ y hiVpityexpl^r J / 

1=1 Pi (t')=p'. 



^j2 h "p"~ AV M{p}) 



1=1 



dt 



Ah 



-I 



dE 
2nK 



Pi{ s )=Pi 



iET/h 



[ ds'U [ M-^V Pl 
Jo ,_i J 



00 



Xexp ^y 



3 -2 



Pi(0)=^ 

ft 2 3 



i^ + es- — Ym 1 r 2 + 2r' 



e i7iT/2m ^ 

47T 2 d COS 1pg(3) 



8, 



2 f-fM t 

^5(3) (q"> q') 



8m 

ftT \ 
2nm J 



i=i 



(is 



(5.112) 
(5.113) 



and ifis(3) (q", q') must be expressed in confocal ellipsoidal coordinates. Actually, "The 
linear element of 5. 2.7.5 and 5.2.7.6 can be considered as degenerations of 5.2.7.7 when 
b = c and a = b, respectively ..." [90] . 



5.3. The Pseudosphere 



5.3.1. General Form of the Propagator and the Green Function. The propagator, 
respectively the Green function, on three-dimensional hyperbolic space, i.e. the three- 
dimensional pseudosphere A^ 3 ) are given by [51, 56, 60] 



z^A (3) ( A i ii i\ rr\ I m \ * rf A( 3)(q // ,q / ) 
K (d A(8) (q ,q);Tj = [ ^-^ j s inhd A(3) (q", qO ^ 



i m 



2RT 4<.»(q",q') 



" 2m" 
(5.114) 



G A(3) (rf A( 3)(q",q');i?) = - 



— m 



2 h 2 sinh<i A (3) (q", q') 

1/2 



x Q 1/2 , _ (coshrf A( 3 ) (q ,/ ,q / )) • (5.115) 

-iy/2mE/h 2 -l-l/2 K AWVM v ; 

where q can denote any of the 34 coordinate systems which separate on A^ 3 ). 



5.3.2. Spherical Coordinates. We consider the coordinate system 

Uo = COsh T\ COsh T2 , Ti , T2 G Hfc , 

■ui = cosh n sinh T2 , 
«2 = sinh Ti sin , < <fi < 2% , 
us = sinh Ti cos (j) . 



i 



(5.116) 



The metric tensor is (g a b) = diag(l, cosh 2 n, sinh 2 n). For the momentum operators we 
have 



ft/ a i 



i 



Pn = y ( -^r + 2 Coth Tl + 2 tanh Tl ' ' 



Pr 2 = , P^ = -qT • (5.117) 

i <3T2 i a0 



The hyperbolic distance is given by 
coshe^^q", q') = coshr{ coshr" cosh(r2 — t' 2 ) — sinhr{ sinhr{ cos(0" — (j)') . (5.118) 
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For the Hamiltonian we have 



h 2 



2m 



M.(s) 



2m 



d 2 



dr 2 \ 



+ ( tanh T\ + coth tA H - 



d 2 1 
+ 



d 



2 1 



>\ P: + coshV/" 2 sinh^n 



) dri cosh 2 n 5r| sinh 2 n <9(/> 2 

+ V 4 + ^ i_ 

8m \ cosh 2 T\ sinh 2 ri 



l 2 \ h 



(5.119) 



We obtain the path integral formulation 

T2{t )=r 2 



Tl(t )=T 1 

/ 

n(t')=ri 



*(t")=*" 



y cosh n sinh nPri(t) ^ £>T 2 (f) ^ V<f>(t) 



T 2 {t')=T^ 



4>(t')=<t>> 



xexp n Jt , 



m 



— ( 7"i + cosh 2 nrf + sinh z ti 
2 



2 _ 12 



) 



ft' 



4 + 



8m V cosh 2 ti sinh 2 ti 



dt 



^ 1/2 

= ( sinh t[ sinh cosh r[ cosh t(') e" i hT/2m ^ 



l^2Z 



2n 



2?r 



x y £>n(t)exp 

n(t')=Ti 





m 9 


-( 


(ft/ 


. 2 Tl " 


2m \ 



,2,1 /2 _ 1 



dt 



= ( sinh t[ sinh r[' cosh r[ cosh t" ] > 



X 



dp*J ,1 ^*(rO*J ,1 ' I) (rOexp 



2tt 
~2m 



JR. 



2tt 



(P 2 + 1) 



(5.120) 



In the path integration, the cylindrical 4>- and free particle T2-path integration have 
been separated straightforwardly. In the last path integral the path integral solutions 
for modified Poschl- Teller potential with v = pi, r\ = \l\ have been used. The wave- 
functions tfjT'Vi) 

are given by 



y&bi> l \ Tl ) = A^ Pl ' Z) (coshTi) ipi - 1 / 2 (sinhri)^- 1 / 2 

X 2^1 



^(|/| + ipi-ip + l),^(|/|+ipi + ip + l);|/| + l;-sinh 2 ri 



jy(pi.O 



1 sinh 7175 / 1 + |/| + ipi + ip\ / 1 + |/| + ipi — ip 



Z I ! V 2tt 5 



(5.121a) 
(5.121b) 
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5.3.3. Horicyclic Coordinates. We consider the coordinate system 
1 fl 



Uq 



2 \y 



+ y + 



ry 2 I ™ 2 

2/ 



Ml = , U 2 = , 

y y 



y>o , 

(x 1 ,x 2 ) = x G IR 2 , 



1 /l 



u z = o\--y 

2 V?/ 



y 



(5.122) 



The metric tensor is (g a b) = H/?/ 2 - For the momentum operators we have p Xl 
— i hd Xl 2 and 

3 



i ^ 2y; • 
The hyperbolic distance is given by 

y' 2 + y" 2 + (x'{ -x[) 2 + (x% - x' 2 ) 2 



(5.123) 



cosh<i A (3)(q // ,q / ) = 
For the Hamiltonian we have 

h 2 . 



2m 



& 2 



2y'y" 



Id d 2 d 2 
+ 



(5.124) 



2m y dy 9x 2 dx\ 

3h 2 

8m 



1 / 

2m 



( 2 i 2 2 i 2 2 \ 

[yp y y + y p Xl +y v X2 ) 



+ 



(5.125) 



The path integral discussion for this coordinate system and its D-dimensional general- 
ization has been discussed in [51]. The path integral is solved by exploiting the path 
integral solution of the radial harmonic oscillator; we cite the result. We obtain the 
identity [51] 



x 2 {t")=x'i 



y(t")=y" x x (t")=x'{ 
y(t')=y' x 1 (t')=x' 1 x 2 (t')=x' 2 



T>x 2 {t) exp 



m x 2 + y 2 3fc 



8m 



dt 



y'y" 



X 



dk 



i ik-(x"-x') 



ir 2 (2tt) 
2 r °° 



71 Jo 



dp p smb. rrpexp 



~2m 



(p 2 + 1) 



K ip (\k\y')K ip (\k\y") 



(5.126) 



5.3.4. Pseudospherical-Conical Coordinates. We consider the coordinate system 
uq = coshr , 

ui = sinh r sn( / u, k) dn(z/, k') 
U2 = sinh r cn(«, k) cn(z/, k') 
us = sinh r dn(«, k) sn(z/, k') 



t > , 
k 2 + k' 2 



(5.127) 
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in the notation of 4.2.3. The hyperbolic distance reads 

cosh d A (3) (q", q') = cosh r" cosh r' — sinh r" sinh r' ^ sn fj," sn fj,' dn i/' dn v' 

+ cn cn \i cn i/' cn v' + dn dn \i sn i/" sn z/ ] . 

(5.128) 

The metric tensor g a b is given by 

(gab) = diag [l, sinh 2 r(k 2 cn 2 fx + /c' 2 cn 2 z/), sinh 2 r(k 2 cn 2 ^ + k' 2 cn 2 z/)] , (5.129) 
and the momentum operators are 



h ( d 

Pt = T V dr + COt T ' ' 



and with p M , as in (4.40). The Hamiltonian has the form 

-2 



(5.130) 



2^ Aa(3) -"2^ 



2 r d 2 B 1 

+ 2cothr— + 



1 



dr 2 



dr sinh 2 r k 2 cn 2 /i + /V 2 cn 2 v \ d/J 2 dv 2 



d 2 d 2 

+ 



1 

2m 



1 



1 



sinh 2 r y/k 2 cn 2 + /V 2 cn 2 v 



(pI+pI) 



a/ /c 2 cn 2 fi + k' 2 cn 2 v 



+ 



1_ 

2m 

(5.131) 



The path integral can be written down yielding 

f.ll\ II /4-"\ I' l-t-H\ I 



J sinh 2 rVr(t) J Vfi(t) J Vv(t){k 2 cn 2 fi + k' 2 cn 2 u) 



T(t>)=T> 

x exp 



A*(*')=M' 



v(t')=v' 



2hj t , 



t 2 + sinh 2 r(k 2 cn 2 fx + k' cn 2 z/)(/i 2 + z> 2 ) 



dt- 



~2m 



oo 



(sinh r' sinh r'O^e-^^^E^^^^K^'O^M 

r(t )=t 



X 



£>r(t) exp 



T(t')=T' 



Z=0 K 

t'' 



m 



h 2 (Z + ±) 



1\2 1 



2m sinh r 



eft 



(sinh r' sinh r'r^^E^^'KJ^TO./.K) 
z=o K 



1 

X — 

7T 



(ippsinh7rp|r(ip + / + 1)| 2 exp 



ihT 
~2m 



(P 2 + 1) 



y<V ip r i/2 (coshr')V ip r i/2 (coshr'') 



(5.132) 
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5.3.5. Hyperbolic-Pseudo-Conical Coordinates. 

uq = coshrnc^, k) nc(v, k') , k 2 + k' 2 = 1 , 



u\ = sinhr , 

U2 = cosh r dc(/i, k) sc(v, k') , 
W3 = coshr sc(/i, /c) dc(^, k') . 



r E IR , 



1 



(5.133) 



Actually, these coordinates encode six different coordinate systems in a similar way as 
discussed in 4.3.5. We omit the details and write down the corresponding path integral 
formulation which reads as follows 



-(t )=- 



v{t")=v" 



J cosh 2 rVr(t) J V/i(t) J Vu(t)(k' 2 ncV + k 2 nc 2 v) 

,(t')=n' v(t')=v' 

t 2 + cosh 2 r(k' 2 nc 2 ^ + k 2 nc 2 v)(fi 2 + z> 2 )l eft - ^ 

J 2m 



T(t')=T' 



ira 



/•OO 

= (coshr' coshr")-^" 1 ^/ 2 -^ / dkAl^')Bl^')A k ^")B k , K {v") 

r(i )=t 

x y X>r(£) e> 



T(t')=T' 



1 



2m cosh 2 r / 

(cosh r' cosh r'T'Yl / ^(/'K^^mM 

e - ihT( P *+i)/2m p -ip ( tanhr /)j»ip_ (tanhr") 

(5.134) 



p sinh 7rp <ip 
2 J]r cosh 2 7r/c + sinh 2 7rp 



m 



rf A( 3 ) (q // ,q / 



2ix\7iT ) sinhrf A( 3 ) (q // ,q / ) 



1m 

2nr 



4(3)(q" 5 q') - 



2m 



(5.135) 



where d A ( 3 ) (q", q') must be expressed in hyperbolic pseudo-conical coordinates. 

5.3.6. Pseudospherical Coordinates. We consider the coordinate system 

-u = cosh r , r > 

tlx = sinh r sin sin <f> , < < 7r , 

it2 = sinh r sin 6> cos , < </> < 27r , 

it3 = sinh r cos 6> . 



(5.136) 



This is the usual three-dimensional (pseudo-) spherical polar coordinates system on 
A( 3 ). We have (g a b) = diag(l, sinh 2 r, sinh 2 r sin 2 6>) and the invariant distance is given 
by 

coshrf A(3 )(q // ,q / ) 

= cosh t" cosh t' - sinh r" sinh r' [ sin 0" sin 0' cos(0" - 0') + cos 0" cos 0'] . (5.137) 
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For the momentum operators we have the operators (4.33,5.130). This gives for the 
Hamiltonian 



2m 



2m 
1 

2m 



8t 2 



, d 

+ 2cothr — + 



i / a 



n 8 

+ cote— + 



i a 2 



dr sinh 2 t\09 2 d$ sin 2 9 dcj) 2 



Pr + ~ 



+ 



+ 



sinh r sinh r sin 9 J 8m V sinh r sinh r sin 



4- 



1 



(5.138) 

In the path integral evaluation one successively uses the path integral solution related 
special case of the modified Poschl- Teller potential. Therefore we obtain the path 
integral formulation (and c.f. [10, 56] for its D-dimensional generalization) 



r(t )=- 



8(t")=6" 

J sinh 2 TVr(t) J sin 9V9(t) J V<j>{t) 

6(t')=6> 



T(t')=T' 

x exp 



4>(t")=4>' 

I 

4>{t')=4>' 



h L, 



m 
~2 



f 2 + sinh 2 t9 2 + sinh 2 r sin 2 



^ 2 ) 



4- 



8m V ' sinh 2 r sinh 2 r sin 2 9 
(sinh r' sinh O^e"^ 2 ™^ ^ Y^tf {6" ,<j>") 



dt 



oo £ 



X 



/,//\ // 

"(* )= r r t // 

/ CT(t)exp Si ( 



2=0 m=-l 
t 



m 



2m sinh r 



T(t')=T' 

(sinhr'sinhr")- 172 ^ ^ y, m *(0', <p')Y l m (9" ,<j>") 

1=0 m=-l 

1 f°° . . . o T i^T 



oo £ 



7V 



x — / dppsinh7rp T(ip + / + I)] exp 



2m 



(P 2 + 1) 



x T^i^ (cosh (cosh r") 



p-l/5 



(5.139) 



5.3.7. Hyperbolic-Spherical Coordinates. We consider the coordinate system 

u = cosh T\ cosh r 2 , T\ e H , 

Ml = sinhri , r 2 > , 

U2 = cosh n sinh r 2 sin , < < 2n , 

«3 = cosh n sinh r 2 cos . 

The hyperbolic distance is given by 

cosh d A (3) (q", q') = cosh r" cosh r{ cosh t 2 ' cosh 

— sinh t'( sinh ri — cosh t[' cosh r{ sinh r 2 sinh T2 cos(0" — 0') . (5.141) 



(0.14UJ 



51 



Here (g a b) = diag(l,cosh ri,cosh r± sinh r 2 ), and for the momentum operators we 
have 

Pn = +tanhTl ) ' ( 5 - 142 ) 

and p T2 ?P0 as m (4.57). For the Hamiltonian we have 



h 



n 1 



2m A<3> 2m 



d 



d 



+ 2tanhn- 1 -5- ■ 

on 9n cosh n V<9t| 



a 2 d id 2 

+ cothr 2 - h 



2m 



2 p; 
Pr, + — r 2 



+ 



+ — 4 + 



cosh n cosh n sinh r 2 / 8m \ cosh n cosh n sinh r 2 



<9r 2 ' sinh 2 r 2 <90 
1 



Therefore we obtain the path integral formulation 



(5.143) 



m* )=ti 



/.ff\ ft 
T 2 (t )=T 2 



4>{t")=4>" 



J cosh 2 nX>ri(t) y sinhT 2 £>T 2 (t) ^ X>0(£) 



n(t')=ri 



T 2 (t')=^ 



<t>(t>)=<t>' 



x exp 



t" r 

/Tfh f o 2 o 2 2 " 2 \ \ 

~2~(ji + cosn r i^2 + cosh n sinh r 2 )J 



8m 



4 + 



1 



cosh T\ cosh ri sinh r 2 



dt 



= ^ cosh r{ cosh t'{) 1 
1 



V x w ^ 2tt 

ie2Z 



x 



1 /" p sinh 7rp dp 

- / 2 2 ex P 

2 7ni cosh 7rpi + sinh 7rp 



dpiPisinh7rpi|r(| + ipi - /) V ipi _i (cosh r 2 )P. pi _i (cosh r 2 ') 



ipi- 



i/iT 
2m 



(P 2 + 1) 



^-,(tanh^)^_ 4 (tanh<) . 

(5.144) 



In the path integral evaluation first the path integral solution of the simple Manning- 
Rosen potential, and second of the simple Rosen-Morse potential was used. 

5.3.8. Equidistant Coordinates. We consider the coordinate system 



uq = cosh n cosh r 2 cosh T3 , ti , r 2 , T3 G IR , 

u\ = sinh n , Ui = cosh n sinh r 2 , 
«3 = cosh ri cosh r 2 sinh T3 . 



(5.145) 



The hyperbolic distance is given by 



cosh d A (3) (q", q') = cosh r{ cosh r[' cosh r 2 cosh r 2 
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+ sinh t% sinh r 2 cos (7-3 — Tg)j — sinh r[ sinh t" . 



(5.146) 



This coordinate system was called in [111] the hyperbolic, respectively the Lobachevskian 
one. Here (g a b) = diag(l, cosh 2 t\, cosh 2 t\ cosh 2 r 2 ), and for the momentum operators 
we have for p Tl as in (5.142), and for p T2 ,Pt 3 as in (4.53). For the Hamiltonian we get 



A(3) 



2m 

1 / 2 
I r 

2m 



2m 







d 



+ 



1 



h2tanhri^ , -j , 

oti dTi cosh n V ° T 2 



d 2 , a 

+ tanh T2 t; 1 - 



c» 2 



1 



+ 



<9t 2 ' cosh 2 t 2 <9t 3 
1 



' J cosh 2 ti cosh 2 ri cosh 2 r 2 / 8m \ cosh 2 t\ cosh 2 t\ cosh 2 r 2 



(5.147) 



Therefore we obtain the path integral formulation 



Tl(t )=T 1 



r 2 (t )=r 2 



T 3(t )=r 3 



cosh T{DTi(t) J cosh T 2 VT 2 (t) J T>T 3 (t) 
Ti(t')=r{ r 2 (t')=r^ r 3 {t>)=T> 3 

{■ r*" r 
— y "2~\ 1 cos ' 1 Tl ^ 2 cos ' 1 Tl cos ' 1 7-2713 J 



4 + 



+ 



8m V cosh 2 ti cosh 2 t\ cosh 2 r 2 



eft 



cosh t[ cosh r" cosh r 2 cosh r 2 



-1/2 



IR 



2^ 



X 



X 



1 

2 Jia cosh 7r/c + sinh 2 npi 
If p sinh 7rp dp 



2 Jir cosh rrpi + sinh 717? 



e - ihT^+i)/2m p - ip i (tanhr ^pip _ i (tanhr^) . 

"" 2 1 (5.148) 



In the path integral evaluation one successively uses the path integral solution related 
special case of the modified Poschl- Teller potential. In Appendix 3 we give the path 
integral solution of the D-dimensional generalization of this coordinate system. 



5.3.9. Hyperbolic-Horicyclic Coordinates. We consider the coordinate system 

cosh r / 1 x 2 \ m 

u = — [- + y + — ) , reJR , 
2 \y y J 

x 

u\ = sinh r , it 2 = — cosh r , x e IR , 

y 

cosh r / 1 x 2 \ 

u 3 = — \ -- y , y > . 

2 \2/ 2/ / 



The hyperbolic distance is given by 



cosh<i A (3)(q // ,q / ) 



y' z + y" 2 + (a" - x') 
2y'y" 



i\2 



(5.149) 



cosh t' cosh t" — sinh r' sinh r" . (5.150) 
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Here (g a b) = diag(l, cosh 2 r/y 2 , cosh 2 r/y 2 ) and for the momentum operators we have 
the operators (4.49,5.142), therefore for the Hamiltonian 



A AO) = 



2m 



2m 
1 
2m 



d 2 n , d y 
+ 2tanhr— + ' 



dr 2 



d d 
+ 



1 



dr cosh 2 T\dx 2 dy 2 
h 2 



cosh r 



{yp 2 y y + y 2 pl) 



+ 



2m 



(5.151) 



We therefore obtain the path integral formulation 

Vy{t) 



r(t")=r" y{t')=y" x(t')=x' 



Vx(t) 



T(t')=T' 

x exp 
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Jo 

r(t")= 

/ 



/ exp 
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h 2 pj+i 
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\TiT 
2m 



I 7 77 poo 

= (cosh r' cosh r")" 1 ^!- / dp e ipo(x "- x,) 

^ JO 

/•oo 

x / dpipisinliTrpiiiripidpolyO-^ipidPoly") 

^0 



x 



p sinh 7rp dp 



2 JiR cosh 7rp2 + sinh np ipi 2 



P-^i (tanhr^Jf (tanhr^) e - iRT (P 2 + 1 )/ 2 "» 



(5.152) 



In this path integral solution the path integral identities from the two-dimensional 
Poincare upper half-plane, and the simple Rosen-Morse potential have been used. 



5.3.10. Horicyclic-Spherical Coordinates. We consider the coordinate system 

„2 ~ 

y >o , 



1/1 r- 
wo = ~ - + y H ) , 

2 v.?/ y 



Ml = 



r sin<^ 
2/ 

1 /l 



r cos (/> 
-u 2 = , r > , 



y 



u 3 = o\ - -y — 

The hyperbolic distance is given by 



< 6< 2tt . 



(5.153) 



UJ , „ „ y' 2 + y" 2 + r ' 2 + r" 2 -2r'r"cos( ( p" - 0') 
coshd A( 3)(q ,q) = — 



(5.154) 
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Here (g a b) = diag(l, l,r 2 )/y 2 and the momentum operators are given by (4.10,5.123). 
For the Hamiltonian we obtain 



2m 



A A (3) = 



h 2 

2^ 



1 



d 2 1 d 2 / d 2 



dy 2 y dy 



Id Id 2 

+ -7T- + 



i 2 i 2 2 i 2/ 2 

= 2m~[ yP y y + y Pr + ^ P <t> 



dr 2 r dr r 2 d<j) 2 
h 2 y 2 3ft 2 



8mr 2 8m 



(5.155) 



This coordinate system is, of course very similar to the horicyclic system and has been 
considered in its .D-dimensional generalization in [51] as well. We just cite the result: 
We obtain the path integral identity 



y(t )=y r(t )=r 

J Vy(t) 



y 



y(t')=y' 
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x 



dppsmhiTpKi p (ky')K[ p (ky") exp 



ihT 
~2m 



(P 2 + 1) 



(5.156) 



In the path integral evaluation after the trivial (p-p&th integration the path integral 
solution of the free radial two-dimensional motion has been used. 



5.3.11. Horicyclic-Spheroidal Coordinates. We consider the coordinate system 
1/1 cosh 2 u — sin 2 u\ 

uo = 77 - + y + , y > o , 

2 \y y J 

cosh u cos v sinh u sin v 
u x = , u 2 = , fx > , 



y 



y 



1/1 cosh fx — sin v 
U3 = 7; - - y 

2 V.?/ y 

The hyperbolic distance is given by 



7T < V < 7T . 



(5.157) 



cosh<i A (3)(q", q') = y' 2 + y" 2 + (cosh 2 fx" - sin 2 v") + (cosh 2 fx' - sin 2 v') 

2y'y" L 

— 2 cosh fx" cosh fx' cos v" cos v' — 2 sinh fi" sinh fx' sin v" sin z/J . 

(5.158) 

Here we have (g a b) = diag(l, sinh 2 ^ + sin 2 z/, sinh 2 // + sin 2 v) together with the 
momentum operator p y (5.123), and for p^ and p v as in (5.45). The Hamiltonian 
has the form 



h 



2m A(3> 2m y 



d 2 1 d 



—z- + 



d 2 d 2 
+ 



dy 2 y dy sinh 2 fx + sin 2 v \ dfx 2 dv 2 
55 



1 

2m 



VPyV + y 2 (sinh 2 fx + sin 2 v) 1 ^ 2 (p 2 l + p 2 )(sinh 2 fx + sin 2 v) 1 / 2 



3ft 2 

8m 



(5.159) 



Using the result form the two-dimensional elliptic coordinate system in we obtain 
the path integral identity 



y(t )=y v(t )=v 
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2m 



(5.161) 



where <i A ( 3 ) (q" q') must be expressed in hyperbolic-spheroidal coordinates. 
5.3.12. Horicyclic-Parabolic Coordinates. We consider the coordinate system 
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The hyperbolic distance is given by 
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(5.162) 



coshrf A(3 )(q // ,q / ) 
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y' 2 + y" 2 + \(ti" 2 + v" 2 ) 2 



+ \(e 2 +v' 2 ) 2 - lie 2 - v ,,2 )(e 2 - v' 2 ) - 2eev''v'] . (s.ies) 



Here we have (g a 6) = diag(l, £ 2 + ?? 2 , £ 2 + ?? 2 ) together with the momentum operators 
(4.22,5.123). The Hamiltonian has the form 
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(5.164) 
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Using the path integral solution corresponding to the two-dimensional parabolic coor- 
dinates in we obtain the identity 

y{t )=y v(t )=v )=£ 

Vy{t) 



y(t')=y" v(t')=V' 

x exp 

ty'y" 



im 
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ihT 
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in the notation of 4.1.5. 



(5.165) 



5.3.13. Confocal Ellipsoidal Coordinates. Let us consider the coordinate system defined 
by 



Ur 



= 



(i = 1,2,3) 



Pi - d Pi - c pi - a Pi-b 
where d<c<ps<b<p2<a<pi (c.f. Section 3.1). Explicitly 

2 (Pi ~d){p2 -d)(p 3 -d) 2 _ (pi - c)(p 2 ~ C)(p 3 - c) 



(5.166) 



Ur 



Ur 



(a - d)(b - d)(c - d) ' 
(pi - a)(p 2 - a)(p 3 - a) 



(a — c)(b — c)(d — c) 
(pi ~ b)(p 2 - b)(p 3 - b) 



(5.167) 



(d-a)(c-a)(b-a) ' 3 "~ {d - b)(c - b){a - b) 
The corresponding line element is given by 
ds 2 1 



dt 2 



(Pi ~ P2)(pl ~ Ps) . 2 , (P2 ~ Ps)(P2 ~ Pi) . 2 , (Pg ~ Pl)(p3 ~ P2) . 2 

Pi ' r> / 1 \ P2 7-) , ; \ P3 



P 4 (pi) ^ P4( P 2) " P 4 (P3) 

(5.168) 

where P 4 (p) = (p — a)(p — b)(p — c)(p — d). Note the similarities as well the differences 
in comparison to the confocal ellipsoidal coordinates on the sphere S^ 3 \ We apply the 
separation formula (3.21). We identify (g a b) = diag(/if, h 2 ,, /13), furthermore 



PM ' (i = 1 ' 2,3) 
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AP A { Pl )P^ P 2)P^P2) 
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P 4 (P1)P4(P3) 



, M 3 = 



Pi ~ P2 



P 4 (P1)P4(P2) 



(5.169) 
(5.170) 

(5.171) 



and obtain the following path integral identity 
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(5.172) 
(5.173) 



and ^(q'W) must be expressed in confocal ellipsoidal coordinates. Actually, the 
relation d < c < . . . is only one possibility to obtain a coordinate system on A^ 3 ). 
We can distinguish eighteen further cases, which fall into two classes with eleven and 
seven coordinate systems, respectively. In the first class emerges when ever two of the 
i e i} = { a 7 ^ c 5 d} are equal. They represent the negative curvature cases in comparison 
to the positive curvature ones as noted in 5.2.6. The second class is obtained by a similar 
consideration as for the pseudo-conical coordinate systems in the two-dimensional case. 
C.f. Section 3.1 for the classification into the four general classes A,B,C,D and [90] for 
the explicit construction. Therefore we obtain a total of 34 different coordinate systems. 
This concludes the discussion of the separable coordinate systems on three-dimensional 
spaces of constant curvature. 



6. Discussion and Summary 

In this paper I have discussed path integrals in two- and three-dimensional spaces of 
constant curvature, i.e. the flat Euclidean spaces 

an d £( 3 ), the spheres and 
S( 3 \ and the pseudospheres A^ 2 ) and A^ 3 \ In many cases I could perform the path 
integration explicitly, thanks to some basic path integral solutions. In the majority of 
the cases, in particular on the pseudospheres A^ 2 ) and A^ 3 \ only an indirect reasoning 
was possible. However, in any case, the explicit expressions of the propagators and 
the Green functions are known in terms of the invariant distance in this space. This 
is possible for any dimension, c.f. [31, 32] for the Euclidean space, c.f. [10, 55] for 
the spheres, and c.f. [10, 56, 107] for the pseudospheres, thus providing coordinate- 
independent expressions. Knowing them gives rise to numerous identities connecting 
the path integral formulations, explicit solutions in terms of the spectral expansions, 
and the coordinate-independent general formulas. 

There are several generalizations possible one can think of. First, one can study 
complex Riemannian manifolds [69, 71] and try to study corresponding path integral 
formulations. At least, the complex Riemannian manifold of constant negative curvature 
has already attracted some attention in the theory of automorphic forms, in order 
to construct a Selberg trace formula [108] similar to the usual one [62, 102, 109]. Here 
also a path integral evaluation is possible [51]. However, it is obvious that complexity 
increases in these more elaborated cases. 

Second, one can introduce potential problems and ask for the separable ones. There 
exist some studies of this kind in the literature, e.g. [49] , where the most important one 
is the Kepler-Coulomb problem in spaces of constant (positive and negative) curvature 
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[9, 47]. However, no systematic study has been done until now. In the case of the 
Coulomb-problem, one is on the one hand interested in the symmetry properties and 
transformations between the bases of the system [17, 79, 95, 97], where the 0(4)- 
symmetry of the hydrogen atom is the best known one and lies at the origin of the 
separability of this problem in four coordinate systems. On the other one is interested 
in the relation [63, 95] of the Kepler-problem and the harmonic oscillator in spaces of 
constant curvature. Whereas in flat space the transformation which relates the Coulomb 
problem in IR 3 and the isotropic harmonic oscillator in IR 4 is the Kustaanheimo-Stiefel 
transformation [23, 24], such a transformation is not known for constant (positive 
and negative) curvature. On and A^ 3 ^ the Coulomb problem and the harmonic 
oscillator are only separable in two coordinate system, the spherical and the sphero- 
conical, respectively, the explict form on E^ 3 \ and A^ 3 ) is given in the following 
table: 





Kepler problem 


Harmonic Oscillator 


£(3) 


a 
r 


m 2 2 
— u) r 
2 


s< 3 > 


a 

- — cot 0i 

ix 


— uj 2 R 2 tan 2 0i 
2 


A (3) 


a 

— — coth r 
R 


VTl o t->2 i 2 

u> R tanh r 

2 



(where we have reintroduced R, and the coordinate notation is as in Section 5). And as 
a third, one can introduce constant magnetic fields on spheres and pseudospheres [44, 
48]. The corresponding cases in flat space are well-known, e.g. [31, 32, 58, 74]. 

A systematic study for the search of separable potentials in IR 2 and IR 3 does 
exist [26, 30], a project which has been started by Smorodinsky, Winternitz et al. 
[81]. These potentials can be put into two classes, which are called maximally and 
minimally superintegrable, respectively. A Hamiltonian system in three degrees of 
freedom is called maximally superintegrable if it admits five globally defined and single- 
valued integrals of motion (the Coulomb-problem and the harmonic oscillator belong 
to this class). A Hamiltonian system in three degrees of freedom is called minimally 
superintegrable if it admits four globally defined and single-valued integrals of motion 
(with the ring-potential as an example [51]). The maximally superintegrable potentials 
can be stated explicitly, whereas the (eight) minimally superintegrable potentials allow 
in addition to an explicit expression an arbitrary function of the coordinates according 
to either oc F(r), oc F(z) or oc F(y/x), respectively. The corresponding super-integrable 
potentials in IR 2 can be obtained in a similar way, and there are a total number of four 
independent potentials [34] . The five maximally superintegrable potentials in IR 3 are: 

V 1 {x,y,z) = -u{x + y +z ) + —^- + - + -j ■ (6.1) 

this potential is separable in (with the coordinates systems which admit explicit path 
integral evaluation in italic) cartesian, spherical, circular cylindrical, elliptic cylindrical, 
oblate spheroidal, prolate spheroidal, conical, and ellipsoidal coordinates. 



V 2 (x,y,z) 



a 



h 2 (k x k 2 



^x 2 + y 2 + z 2 2m\x 2 y 2 



(6.2) 
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this potential is separable in spherical, parabolic, and conical coordinates. 



K,(*,*z) = f (-^ + *| + *|) ; (6.3) 
this potential is separable in spherical and parabolic cylindrical coordinates. 

^,^ = r ^= - + ; 6.4 

this potential is separable in parabolic cylindrical and parabolic coordinates. 



F 5 (x, y, z) = yW V + J + iz>) + — ^ + -±j ; (6.5) 

this potential is separable in cartesian, parabolic (one obtains an intractable quadratic 
plus sextic oscillator), and elliptic cylindrical coordinates. 

In each case it is easy to construct the various path integral formulations, based 
on the coordinate systems as listed in Section 5.1, and obtain therefore the various 
identities connecting them. The corresponding path integral evaluations then can be 
performed by using the basic path integrals as given in Appendix 1, which is left as an 
exercise to the interested reader. Note that in comparison to the free motion in spaces 
of constant curvature the statement of the propagator in a coordinate independent way 
is not possible. This concludes the discussion. 
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Appendix 1. Some Important Path Integral Solutions 

In this Appendix we cite some important path integral solutions, in particular for the 
radial harmonic oscillator and for the (modified) Poschl- Teller potential, including some 
special cases. These three path integral solutions are the most important building blocks 
for almost all other path integrals. 

A. 1.1. The path integral for the radial harmonic oscillator. 

The calculation of the path integral for the radial harmonic oscillator has first been 
performed by Peak and Inomata [92]. A more general case is due to Goovaerts [37] 
(c.f. also [22]). Path integrals related to the radial harmonic oscillator may be called of 
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Besselian type [65] . Here we are not going into the subtleties of the Besselian functional 
measure due to the Bessel functions which appear in the lattice approach [33, 55, 104] 
which is actually necessary for the explicit evaluation of the radial harmonic oscillator 
path integral [21, 37, 92]. One obtains (modulo the mentioned subtleties) (r > 0) 
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(A1.1) 



In the special case that a; = 0, i.e. for the free radial motion, we obtain 

r{t")=r'' 



(A1.3) 
(A1.4) 



A. 1.2. The modified Poschl-Teller Potential. 

The path integral solution for the Poschl-Teller potential can be achieved by means of 
the SU(2)-path integral. We have [10, 21, 66] 
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with the wavefunctions given by 
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A.1.3. The modified Poschl-Teller Potential. 

The path integral solution for the modified Poschl-Teller potential can be achieved by 
means of the SU(1, l)-path integral. We have [10, 66, 68] 
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J Vr(t) exp 
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+ ^"dp^)*(rO^ ) (Oexp^-^^ . (A1.7) 

Now introduce the numbers fci, fc 2 defined by: k\ = |(1 ± u), k 2 = |(1 ± 77), where the 
correct sign depends on the boundary-conditions for r — > and r — > 00, respectively. In 
particular for rf = ^, i.e. fc 2 = j, |, we obtain wavefunctions with even and odd 
parity, respectively. The number A/m denotes the maximal number of states with 
0, 1, ... , N M < ki — k 2 — \. The bound state wavefunctions read as (k = k\ — k 2 — n) 



$(*i,*2)( r ) = ^ fcl ' fc2 )(sinhr) 2fc2 -i(coshr)- 2fcl+ i 
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(A1.8) 



The scattering states are given by: 



$(fci,fc 2 )( r ) = ^ fcl ' fc2 )(coshr) 2fcl -5(sinhrj 
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> (A1.9) 



x r(fei + k 2 + K- i)r(-fei + k 2 - k + 1) 



[k = |(1 + ip)]- Of course, in the path integral formulation of the modified Poschl-Teller 
potential a similar functional weight interpretation must be used as for the Poschl- 
Teller potential in order to have a proper short-time behaviour, respectively a lattice 
regularization [46]. 

I also cite two special cases, where first only the 1/ sinh 2 r, and second where only 
the 1/ coshr term is present. The special case V^ sh \r) = Vq/ sinh 2 r (Vo = ^(A 2 — \)i 
r > 0) is given by [46, 56, 75] (simple Manning-Rosen potential) 
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= — a/suiIi r' sinh r" / <ippsinh7rp|r(i + ip — A)| 2 
^ JO 

x Irll + ip-^l^^^coshrOPr-i^^oshr'Oe- 1 ^ 2 ^ 2 - . 

(ALIO) 

The special case V^ ch \x) =V / cosh 2 x (V = |^(& 2 + \), x G IR) on the other is given 
by [45, 46, 75] (simple Rosen-Morse potential) 
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Note that for Vo < also bound state solutions are allowed. These, however, we do not 
need. 



Appendix 2. Discussion of a Dispersion Relation 

We consider the integral representation as given by Buchholz [13, p. 158]: 
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(A2.1) 



In order to apply (A2.1) for the determination of continuous spectra we perform 
some manipulations. We replace x — > 2x, y — > — 2y, s — > + ip, Xi,2 — ¥ (1 + /-0/2, 
sin</> — > 1/ sin a, and — ► This gives 
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(A2.2) 



where use has been made of some properties of the Whittaker functions M^ v {z) = 
r(l + 2/i)M M ,„ and [13, p.ll] M x ^ /2 (ze ±i ^) = e ±i7r ( 1+ ^/ 2 M_ Xi/x/2 (z), respectively. 
Equation (A2.2) is the desired expansion formula for the determination of continuous 
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spectra. Note the additional e^-factor in (A2.2). This representation can also be 
deduced from the integral formula [29, p. 414; 41, p. 884] 
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We consider the special case of (A2.2), where the Bessel function is replaced by an 
exponential, i.e.: 
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Together with the properties of the Whittaker functions that for ji = ±| they can be 
expressed by parabolic cylinder functions eI 0,1 ^ 1 we obtain 
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The functions £^ 0) (» and E^ (z) are even and odd in the variable z, respectively. 
We can establish the connection to the parabolic cylinder functions D v and (A2.4) is 
equivalent with the expansion [41, p. 896]: 
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Appendix 3. The D-Dimensional Hyperbolic System 



We consider the D-dimensional generalization of the hyperbolic coordinate system: 



Uq = cosh n .. .cosh Td_i , 

u\ = sinhri , 

U2 = cosh 7~i sinh r 2 , 

ud-i = coshri . . . sinhrD_i . 



(A3.1) 



Therefore we obtain the path integral formulation by successively applying (Al.ll) 
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(A3.2) 



with the wave-functions ^^(t) given by 



PsinhTrp pi v t tanh _\ 

2 (cosh 2 ttA; + sinh 2 np) lfc 2 



(A3.3) 



This solution represents an alternative path integral solutions as already outlined for 
two other coordinate systems in [51]. Note that irrespective which coordinate system 
for the path integral formulation on A^ -1 ) is chosen (and there are many indeed) the 
Green function always has the form [51, 56] 
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